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wise.  The  periodic  time-dependent  nature  of  the  bunching  results  in  effi¬ 
cient  microwave  generation.  Experimental  efficiencies  of  12%  have  been 
obtained  in  previous  work  and  theoretical  efficiencies  of  20%  have  been 
obtained  in  the  work  described  here.  The  microwave  frequency  scales  as  the 
relativistic  beam  plasma  frequency.  This  can  be  adjusted  via  an  externally 
applied  axial  magnetic  field.  The  result  is  an  oscillator  which  is  fre¬ 
quency  agile  over  an  order  of  magnitude  without  externally  changing  the 
physical  configuration.  Output  powers  of  from  kilowatts  to  gigawatts  are 
obtainable  based  on  the  power  source.  Experimental  results  are  reviewed 
and  explained.  The  rationale  for  using  a  foil  less  diode  as  the  electron 
beam  source  is  given. 

A  detailed  mechanical  and  electrical  design  has  been  presented  in  pre¬ 
vious  reports  (AMRC-R-451,  1983).  It  entails  placing  the  vi rcator  in  a 
high  vacuum  vessel  with  flexibility  for  further  development.  The  electri¬ 
cal  power  is  pulse  formed  using  a  lumped  artificial  Blumlein  line.  The 
axial  magnetic  field  is  critically  damped  for  component  safety  and  circuit 
design  simplicity.  The  beam  source  is  a  foilless  diode  with  an  annular 
limiter  for  collimation  and  prevention  of  electrons  reflexing  into  the 
diode. 

Several  important  technical  developments  have  come  to  a  successful 
conclusion  during  the  period  of  work  covered  by  this  report.  The  multi- 
Tesla  magnetic  field  coil,  critically  important  for  achieving  the  high 
electron  densities  necessary  for  driving  millimeter  wavelenghts,  has  been 
designed,  and  successfully  tested.  Secondly,  a  millimeter  wavelength  spec¬ 
trometer,  necessary  for  diagnosing  the  time  history  and  spectral  content  of 
the  vi rcator  emission,  has  been  designed  and  calibrated  over  a  wide  band 
from  35  to  110  GHz.  Measured  frequency  selectivity  was  typically  of  order 
5%  or  better,  i.e.,  af/f  <  .05.  A  unique  cathode  material  has  been 
investigated,  yielding  fast  field-emission  turn  on  and  low  gap  closure 
rates.  Details  of  these  results  are  presented  in  the  main  section  of  this 
report. 
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SECTION  I 
INTRODUCTION 


The  main  thrust  of  research  in  the  second  year  of  this  program  has 
been  experimental.  This  report  summarizes  the  progress  made  in  developing 
the  hardware  for  a  millimeter  wavelength  virtual  cathode  microwave  oscil¬ 
lator  (Vircator).  Fabrication,  testing,  and  calibration  of  two  critical 
components,  the  solenoidal  magnet  and  the  millimeter  microwave  spectrom¬ 
eter,  is  complete.  Data  summary  and  complete  documentation  of  the  testing 
of  the  magnetic  field  coil  at  the  30  kG  level,  is  presented  in  Section  II. 
Calibration  data  for  the  spectrometer,  covering  the  entire  35  to  100  GHz 
frequency  range,  is  presented  in  Section  III.  Section  IV  describes  the 
successful  operation  of  the  electron  beam  pulse  power  system,  based  on  an 
artificial  Blumlein  design.  A  unique  cathode  material  is  described  which 
allows  fast  turn  on  of  electron  field  emission,  at  field  strengths  of 
approximately  50  kV/cm  and  very  slow  anode-cathode  gap  closure  rates,  less 
than  1  cm/usec.  Our  work  in  this  area  is  further  supplemented  by  reprints 
of  articles  compiled  in  Appendix  B. 


SECTION  II 
MAGNET  DESIGN 


The  main  criteria  for  the  vircator  magnet  design  are: 

(1)  Peak  magnetic  field  of  60  kG; 

(2)  Field  homogeneity  in  the  anode-cathode  gap  region  of  aB/B  <  1%. 

(3)  Accessibility  to  the  high  magnetic  field  region  for  grids,  elec¬ 
trodes,  diagnostics,  and  their  connecting  cables  and  mechanical 
support  structure. 

In  support  of  the  vi rcator  magnet  design,  a  useful  computational  tool 
has  been  developed  for  calculating  the  field  line  configuration  for 
arbitrary  solenoidal  geometries.  The  computer  code  uses  a  two  step 
process.  Radial  and  axial  components  of  the  magnetic  field  are  calculated 
on  a  two  dimensional  rectangular  grid  by  summing  the  weighted  contribution 
from  each  circular  current  loop  making  up  the  solenoid.  Bilinear 
interpolation  is  then  used  to  follow  field  line  trajectories  which  are 
launched  from  the  center  of  the  solenoid. 

The  magnet, c  field  for  an  arbitrary  solenoidal  configuration  follows 
from  the  weighted  superposition  of  the  fields  from  the  elemental  circular 
current  loops,  with  the  neglect  of  helicity  and  finite-current  cross- 
section  effects.  Eddy  current  effects  are  also  neglected;  this  is  a  good 
assumption  since  the  skin  penetration  time  through  the  high  resistance 
stainless  steel  casing,  approximately  10  microseconds,  is  small  relative  to 
the  570  microsecond  field  rise-time.  With  these  approximations,  the  radial 
and  axial  components  of  the  magnetic  field,  Br  and  Bz,  at  an  arbitrary 
point  r,z,  due  to  a  circular  current  loop  of  radius  a  centered  at  the 
origin  of  a  cylindrical  coordinate  system  (see  Figure  1)  can  be  calculated 
(in  mks  units)  as1 : 
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Figure  1.  Relation  of  the  cylindrical  co-ordinate  system  to  the 
circular  current  loop.  The  center  of  the  loop  and  the 
origin  of  co-ordinates  are  coincident,  and  the  z  axis 
and  the  normal  to  the  loop  are  coincident. 
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and  K(k)  and  E(k)  are  the  complete  elliptic  integrals  of  the  first  and 
second  kind,  respectively.  It  is  straightforward  to  evaluate  the  elliptic 
integrals  numerically,  and  sum  over  the  windings  of  the  solenoid  to  produce 
a  Br  and  Bz  matrix  representing  the  calculated  magnetic  field.  The 
first  part  of  the  field  line  plotting  carries  out  these  calculations,  and 
is  based  on  a  code  provided  by  John  Freeman  of  Sandia  National 
Laboratories.2  Run  time  is  approximately  30  seconds  on  a  VAX  11/780  for  an 
84  point  (12x7)  grid  and  a  156  turn  solenoid. 


Field  line  plots  are  produced  from  the  R  field  matrix  using  a  shooting 
technique.  Field  lines  are  launched  from  the  center  of  the  solenoid  in  the 
x-y  plane,  nominally  a  uniform  field  region,  and  the  trajectory  incremented 
by  step  sizes  which  are  small  relative  to  the  B  field  grid  points. 

Bilinear  interpolation  allows  calculation  of  Br  and  Bz  between  grid 
points  with  second  order  accuracy.3 


In  order  to  verify  accurate  operation  of  the  code,  a  cross  check  has 
been  carried  out  with  an  analytical  calculation  of  a  field  line  trajectory 
near  the  magnet  center.  The  magnetic  field  in  the  central  zone  of  a  uni¬ 
form  current  density,  cyl indrical ly  symmetric  solenoid  can  be  written  as  a 
power  series  with  terms  of  the  form 
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where  J0,  the  average  solenoid  current  density,  has  units  of  amps/cm^, 
magnet  radius  a  is  in  cm,  and  B0  is  in  gauss.  The  equation  defining  the 
field  line  trajectory  in  two-dimensional  space  can  be  written: 

dr  _  dz 


Keeping  only  the  first  term  in  the  B  field  expansion  and  integrating  the 
resulting  differential  equation  gives  the  equation  for  the  field  line  tra¬ 
jectories  near  the  center  of  solenoid: 


r  =  rQ  e 


■E2z2/2a2 


For  the  case  of  a  10  cm  radius,  100  cm  long  solenoid  wound  with  a  sin¬ 
gle  251  turn  layer  (a  *  1,  B  =  5,  E2  *  -0.0024),  the  code  calculation  shows 
that  a  field  line  launched  from  (r,z)  =  (0.90000,  0)  cm  intersects  a  point 
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(0.90095,  9.375)  cm.  Taking  r0  =  0.9  in  Equation  (17)  yields  (r,z)  = 
(0.9009497,  9.375  cm),  which  agrees  with  the  numerical  calculation  to 
within  a  fraction  of  one  percent.  In  addition,  the  B  field  intensity 
at  the  center  of  the  magnet  for  1  anp  current  is  calculated  numerically  as 
3.081  Gauss,  in  excellent  agreement  with  an  analytical  calculation8 
yielding  3.0806  G  for  this  geometry.  These  cross-checks  indicate  that 
both  the  IT  field  and  the  field  line  trajectory  calculations  are  working 
correctly  for  the  case  of  interest,  i.e.,  near  the  center  of  a  long 
solenoid. 

The  magnetic  field  code  has  been  the  primary  tool  used  in  design  of 
the  vircator  magnet.  End  windings  were  used  to  yield  a  homogeneous  central 
field.  Optimization  was  carried  out  by  iterating  on  a  base  design  using 
the  interactive  capabilities  of  the  VAX  11/780.  A  listing  of  the  FORTRAN 
source  code  is  given  in  Appendix  A.  Table  I  summarizes  the  final  magnet 
design  parameters,  and  Table  II  gives  the  radial  and  axial  magnetic  field 
components  over  a  central  region  3  cm  in  diameter  and  8  cm  long.  Homogen¬ 
eity  in  this  region  is  aB/B0  _<  0.008,  Br/B0  <_  0.007,  with  B0  the 
field  at  the  center  of  the  magnet. 

The  magnet  overall  length  is  21.4  cm  with  a  usable  inside  diameter  of 
4.6  cm.  The  magnet  is  wound  with  156  turns  of  1  mm  diameter  wire  in  5 
layers,  yielding  an  inductance  of  0.7  mH  and  a  current  of  5.7  kA  for  the 
full  60  kG  field.  Self-forces  tending  to  squeeze  the  magnet  axially  are 
taken  up  by  the  copper  windings  and  an  epoxy  binder.  Radial  expansion 
forces  are  taken  up  by  an  encapsulating  cylindrical  stainless  steel  tube 
and  the  epoxy  binder.  With  the  assumption  that  half  the  radial  expansion 
forces  are  borne  by  the  epoxy  and  half  by  the  encapsulating  stainless 
steel,  the  stress  in  the  1.9  mm  thick,  11.4  cm  O.D.  stainless  shell  is 
23,100  psi  for  a  60  kG  field.  Since  the  yield  strength  of  type  304  stain¬ 
less  steel  exceeds  40,000  psi,  the  design  provides  an  acceptable  safety 
factor. 


15 


TABLE  I.  VIRCATOR  MAGNET  CONFIGURATION 


Shell 

Number 

Conductor 

Radius 

(cm) 

Conductor  Axial 
Coordinate  Range 
(cm) 

Number  of  Windings 
Distributed  in 
Shell 

Description 

1 

2.80 

-7.41  to  -4.94 

7 

Inner  end-correc¬ 
tion  winding, 
left  side. 

2 

2.80 

+4.94  to  +7.41 

7 

Inner  end-correc¬ 
tion  winding, 
right  side. 

3 

3.27 

-7.41  to  -0.41 

17 

Outer  end-correc¬ 
tion,  left  side. 

4 

3.27 

+0.41  to  +7.41 

17 

Outer  end-correc¬ 
tion,  right  side. 

5 

3.73 

-7.41  to  +7.41 

36 

Main  winding. 

6 

4.20 

-7.41  to  +7.41 

36 

Main  winding. 

7 

4.66 

-7.41  to  +7.41 

36 

Main  winding. 

Design  Magnetic  Field  .  60  kG 

Total  Number  of  Turns  . . .  156 

Drive  Current  for  60  kG  Field  at  Center  . . .  5.7  kA 

Magnet  Inductance  .  0.7  mH 

Field  Energy  at  Full  Current  . . .  11.4  kJ 


TABLE  II.  MAGNET  COMPUTATION  REGION  DATA 


RADIUS=  1.500E+00 
DZ=  l.OOOE+OO 


LENGTH=  4.000E+00  DR3  5.000E-01 
NR=  4  NZ=  5 


COIL  DATA 
R=  2.80E+00 
R=  2.80E+00 
R=  3.27E+00 
R=  3.27E+00 
R3  3.73E+00 
R=  4.20E+00 
R=  4.66E+00 


ZL=-7 .41E+00 
ZL=  4.94E+00 
ZL=-7 .41E+00 
ZL=  4.12E-01 
ZL=-7 .41E+00 
ZL=-7.41E+00 
ZL=-7 .41E+00 


ZR=-4.94E+00 
ZR=  7.41E+00 
ZR=-4.12E-01 
ZR=  7.41E+00 
ZR=  7.41E+00 
ZR=  7.41E+00 
ZR=  7.41E+00 


TURNS3  7.00E+00 
TURNS3  7.00E+0u 
TURNS3  1.70E+01 
TURNS3  1.70E+01 
TURNS3  3.60E+01 
TURNS3  3.60E+01 
TURNS3  3.60E+01 


CURRENT3  5.71E+03 
CURRENT3  5.71E+03 
CURRENT3  5.71E+03 
CURRENT3  5.71E+03 
CURRENT3  5.71E+03 
CURRENT3  5.71E+03 
CURRENT3  5.71E+03 


RADIAL  MESH  LINE  LOCATIONS 

O.OOOOOE+OO  5.00000E-01  l.OOOOOE+OO  1.50000E+00 


AXIAL  MESH  LINE  LOCATIONS 

O.OOOOOE+OO  l.OOOOOE+OO  2.00000E+00  3.00000E+00  4.00000E+00 

THETA  =  ARCTAN  (BR/BZ)  (degrees) 

LISTING  OF  BR  AND  BZ  OVER  THE  COMPUTATIONAL  GRID 


K 

L 

R 

z 

BR(K.L) 

bz(k,l) 

B(K,L) 

BR/BZ 

THETA 

1 

1 

0.00000 

0.0000 

0.000 

60000.031 

60000.031 

O.OOOE+OO 

0.0000 

1 

2 

0.00000 

1.0000 

0.000 

60017.418 

60017.418 

O.OOOE+OO 

0.0000 

1 

3 

0.00000 

2.0000 

0.000 

59966.418 

59966.418 

O.OOOE+OO 

0.0000 

1 

4 

0.00000 

3.0000 

0.000 

59726.496 

59726.496 

O.OOOE+OO 

0.0000 

1 

5 

0.00000 

4.0000 

0.000 

59156.402 

59156.402 

O.OOOE+OO 

0.0000 

K 

L 

R 

Z 

BR(K,L) 

BZ(K.L) 

B(K,L) 

BR/BZ 

THETA 

2 

1 

0.50000 

0.0000 

0.000 

59995.840 

59995.840 

6.835E-10 

0.0000 

2 

2 

0.50000 

1.0000 

-4.416 

60021.887 

60021.887 

-7.357E-05 

-0.0042 

2 

3 

0.50000 

2.0000 

33.265 

59978.754 

59978.762 

5.546E-04 

0.0318 

2 

4 

0.50000 

3.0000 

87.598 

59743.328 

59743.391 

1.466E-03 

0.0840 

2 

5 

0.50000 

4.0000 

214.361 

59216.949 

59217.336 

3.620E-03 

0.2074 

K 

L 

R 

Z 

BR(K,L) 

BZ(K.L) 

B(K,L) 

BR/BZ 

THETA 

3 

1 

1.00000 

0.0000 

0.001 

59978.680 

59978.680 

1.345E-08 

0.0000 

3 

2 

1.00000 

1.0000 

-18.919 

60036.184 

60036.184 

-3.151 E-04 

-0.0181 

3 

3 

1.00000 

2.0000 

65.006 

60017.195 

60017.230 

1.083E-03 

0.0621 

3 

4 

1.00000 

3.0000 

166.902 

59787.082 

59787.316 

2.792E-03 

0.1599 

3 

5 

1.00000 

4.0000 

365.696 

59381.059 

59382.184 

6.158E-03 

0.3528 

K 

L 

R 

z 

BR(K,L) 

BZ(K,L) 

8(K,L) 

BR/BZ 

THETA 

4 

1 

1.50000 

0.0000 

0.000 

59929.898 

59929.898 

4 .424E-09 

0.0000 

4 

2 

1.50000 

1.0000 

-58.797 

60064.949 

60064.977 

-9.789E-04 

-0.0561 

4 

3 

1.50000 

2.0000 

96.956 

60087.793 

60087.871 

1.614E-03 

0.0925 

4 

4 

1.50000 

3.0000 

241.009 

59839.746 

59840.230 

4.028E-03 

0.2308 

4 

5 

1.50000 

4.0000 

391.128 

59583.676 

59584.961 

6.564E-03 

0.3761 

m 


K  * 
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Encapsulation  of  the  magnet  is  completed  by  a  1.6  mm  thick  inner  wall 
and  two  1.6  mm  thick  end  plates.  The  stainless  shell  and  its  relation  to 
the  coil  windings  is  shown  in  Figure  2.  The  inner  stainless  cylinder  acts 
to  shield  the  magnet  winding  from  electron  and  x-ray  bombardment  which 
would  otherwise  reduce  the  magnet  turn-to-turn  voltage  standoff.  Since  the 
highest  operating  voltage  for  the  vircator  is  100  kV,  the  inner  stainless 
shield  provides  complete  shielding  from  electron  bombardment  (50  electron 
ranges  at  100  kV)  and  better  than  a  factor  of  10  attenuation  (2.4  mean  free 
paths  for  50  kV  x-rays)  for  the  most  energetic  x-rays  of  interest.  Flux 
penetration  constraints  do  not  allow  the  wall  thickness  to  be  substantially 
increased  beyond  2  mm,  corresponding  to  a  magnetic  field  penetration  time 
of  approximate! j  10  microseconds. 

Fabrication  prints  for  the  stainless  steel  magnet  housing  are  shown  in 
Figures  3-5.  Fabrication  prints  for  the  PVC  insulators  are  shown  in 
Figures  6-9,  and  the  assembled  magnet  is  shown  in  Photo  1.  The  purpose  of 
the  insulator  assembly  is  to  provide  high  voltage  standoff  between  the 
multi-kV  magnet  windings  and  the  grounded  stainless  steel  magnet  housing. 
PVC  wall  thickness  is  5.7  mm  on  the  outer  wall  and  24.9  mm  on  the  side 
walls,  providing  voltage  standoff  to  well  above  the  10  kV  required.  No  PVC 
is  required  at  the  inner  radius  since  the  inner-most  magnetic  winding  is 
relatively  close  to  ground  potential. 

The  cylindrical  magnet  is  held  in  place  by  encapsulation  between  the 
upper  and  lower  PVC  brackets  shown  in  Figure  10.  Brackets  attach  to  the 
upper  plate  of  Figure  11  which  is  positioned  with  respect  to  the  lower 
plate  of  Figures  12  and  13  by  3  bolts.  These  bolts  provide  3  degrees  of 
freedom  for  aligning  the  magnet  with  the  cathode.  An  additional  two 
degrees  of  freedom  is  allowed  by  positioning  the  lower  plate  of  Figure  12 
relative  to  the  stainless  steel  mounting  brackets  which  are  welded  inside 
the  vacuum  vessel;  both  of  these  pieces  have  racetrack-type  slots  for 
adjustment.  Thus,  full  5  degree-of-f reedom  flexibility  is  allowed  by  the 
magnet  alignment  system.  Photo  2  shows  the  vircator  magnet  and  the 
assembled  alignment  system. 


t 
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Figure  3.  Fabrication  print  of  the  magnet  outer  shell  used  to  take 

up  the  radial  stresses  induced  by  the  60  kG  magnetic  field. 


ation  print  of  the  PVC  spacer  separating 
correction  windings. 
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Photograph  1 . 


Vircator  magnet  assembly.  The  cylindrical  magnet  housing 
is  attached  to  the  PVC  baseplate  with  four  stainless  steel  rods. 
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UPPER  PLATE 


Fabrication  print  of  the  PVC  upper  plate  which  provides  a  rigid  mount 
for  the  magnet. 


Figure  13.  Fabrication  print  of  the  PVC  pieces  used  to 
stiffen  the  lower  plate  structure. 


Photograph  of  the  vircator  magnet  and  the  five  degree-of- 
t'-eedo”  positioning  assembly  for  aligning  the  magnet  with 
tnn  high  voltage  cathode. 
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A  schematic  for  the  circuit  which  drives  the  magnetic  field  coil  is 
shown  in  Figure  14.  A  366  microfarad  energy  storage  capacitor  bank  is  dis¬ 
charged  into  the  magnetic  field  coil  through  a  series  resistance  which  is 
chosen  so  that  the  circuit  is  strongly,  but  not  critically,  damped.  Since 
the  current  reversal  is  slight,  the  lifetime  of  both  the  capacitors  and  the 
magnetic  field  coil  are  increased.  The  price  paid  is  an  increase  by  a  fac¬ 
tor  of  approximately  2  in  the  required  bank  capacitance  relative  to  an 
undamped  circuit.  The  same  result  could  be  achieved  by  crowbaring  the  cur¬ 
rent  at  peak  field  with  an  ignitron,  and  eliminating  the  series  resis¬ 
tance.  Since  the  cost  of  the  increased  capacitance  is  small  compared  to 
the  cost  of  the  ignitron  and  its  driving  and  control  circuits,  the  former 
option  was  chosen. 

The  magnet  energizing  circuit  of  Figure  14  operates  as  follows.  A 
20  kV,  10  mA  power  supply  (Del  Electronics  IRIU-20-10)  charges  the  366 
microfarad  capacitor  bank  to  its  charge  voltage,  about  8  kV  for  a 
30  kilogauss  field.  A  krytron  gas  tube  switch  is  command  triggered  to  dis¬ 
charge  a  1  microfarad  capacitor  charged  to  2.5  kV  into  the  ignitor 
terminals  of  a  mercury  ignitron  {General  Electric  #GL-7703).  The  triggered 
ignitron  then  discharges  the  main  capacitor  bank  into  the  inductive  load 
through  a  1.8  ohm  series  resistor.  The  peak  magnetic  field  occurs  at  570 
microseconds.  A  time  delay  generator  (California  Avionics  #113CR)  command 
fires  the  spark  gap  to  initiate  the  electron  beam  at  peak  magnetic  field 
time. 


Testing  of  the  vircator  magnet  has  been  carried  out  to  the  30  kG  level 
using  two  183  pF,  8  kV  rated  capacitors  configured  in  parallel  as  the 
energy  source.  No  attempt  has  been  made  to  exceed  this  field  level, 
although  design  calculations  indicate  reliable  operation  to  60  kG.  Time 
dependence  of  the  magnetic  field  has  been  measured  with  the  probe  assembly 
shown  in  Photo  3.  A  100  turn  coil,  inside  diameter  9.538  mm,  outside 
diameter  11.15  mm,  and  overall  length  3.17  mm,  is  wound  on  a  PVC  rod  to 
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Fiqure  14  Schematic  of  the  pulse  power  circuit  which  drives  the  vircator  magnet  coil. 

The  series  1.8  resistor  strongly  damps  the  magnet  RLC  circuit  to  minimize 

current  reversal. 


100  turn  coil  is  at  the  far  right.  Measured  values  are  consisten 
with  the  elementary  RLC  circuit  theory  described  in  the  text. 

The  probe  axial  stroke  is  6  inches. 


sample  the  Bz  field.  Because  of  the  slow  field  rise  time  and  low  probe 
inductance,  the  probe  is  not  self  integrating  and  produces  a  Bz  signal. 

The  peak  Bz  field  follows  from  the  stored  voltage  output  versus  time  V(t) 
according  to: 

c  t 

B  (t)  =  1.19  x  10°  /  dt'V(f)  .  (18) 

o 

The  probe  can  be  positioned  anywhere  along  the  axis  of  the  magnet  using  the 
vacuum  manipulator  of  Photo  3  (Huntington  Mechanical  Laboratories  linear 
motion  feedthrough  #VF-158).  A  typical  probe  output  oscilloscope  trace  is 
shown  in  Figure  15  for  a  4.0  kV  capacitor  charge  voltage  with  the  probe 
located  on  axis.  The  magnetic  field  peaks  at  570  psec,  and  carrying  out 
the  areal  integral  in  the  above  equation  yields  a  peak  field  of  14.62  kG. 
The  measured  field  level  and  risetime  can  be  compared  with  analytical  cal¬ 
culations9  based  on  a  simple  series  RLC  circuit.  These  calculations  show 
that  for  a  capacitor  C  initially  charged  to  a  voltage  V0  and  discharged 
at  time  t=0  into  an  underdamped  series  load  of  inductance  L  and  resistance 
R  (mks  units),  the  current  peaks  at  time  t0  at  a  value  Imax  given  by: 


t  =  -  tan'1  (  )  , 

o  to  R  ; 

(19) 

V 

I  _  ° 

max  /L/C 

(20) 

2  1-d 

oj  =  lc 

(21) 

R2C 

d  5  4L~  • 

(22) 

and 
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(23) 


Substituting  the  system  impedances  for  the  vircator  magnet  (R  =  1.8  a,  L  = 
0.7  mH,  C  =  366  yF,  charge  voltage  V0  =  4  kV)  yields  t0  =  575  ysec  and 
Imax  =  1.381  kA.  The  series  resistor  in  this  circuit  is  a  copper  sulfate 
based  liquid  type  with  the  capacity  to  dissipate  20  kj  with  less  than  a 
20°C  temperature  increase.  The  resistance  value  is  determined  by  a  piston 
type  plunger  and  was  not  accurately  measured;  the  1.8  n  value  used  here  was 
chosen  to  give  good  agreement  with  experiment.  The  magnetic  field  at  the 
center  of  the  magnet  B0  is  linearly  proportional  to  the  series  current 
I0  according  to  (Table  II  shows  B0  =  60  kG  for  I0  =  5.71  kA): 

Bq  =  10.51  I 0( kA )  kG.  (24) 

Thus,  for  4  kV  charge  voltage,  the  calculated  field  is  14.51  kG,  and  the 
risetime  575  ysec.  These  values  are  in  good  agreement  with  the  measured 
values  of  14.62  kG  and  570  usee,  indicating  the  magnet  operation  is  well 
characterized  by  a  strongly  damped  series  RLC  circuit. 
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SECTION  III 

CALIBRATION  OF  THE  35-110  GHz  VIRCATOR  SPECTROMETER 
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A  useful  device  for  characterizing  the  intensity  and  wavelength  of 
electromagnetic  radiation  is  the  spectrometer.  Although  Fabry-Perot  inter¬ 
ferometry  and  Fourier  transform  spectroscopy  are  commonly  used  in  the  UV, 
visible,  and  IR  regions  of  the  spectrum,  diffraction  grating  spectroscopy 
carries  over  conveniently  to  microwave  frequencies.  The  fast  response  of 
crystal  detectors  which  can  easily  be  coupled  to  a  microwave  grating  spec¬ 
trometer  makes  this  a  good  choice  for  monitoring  radiation  from  pulsed 
sources.  A  grating  spectrometer  specifically  designed  to  monitor  milli¬ 
meter  waves  has  been  fabricated.  The  design  philosophy  has  been  previously 
described.10’11  Spectrometer  calibration  is  described  here. 

The  main  tool  used  in  calibrating  the  vircator  spectrometer  is  a 
Micro-Now  Model  #7058  millimeter  wave  sweeper  system  with  appropriate  heads 
for  B-band  (35-50  GHz),  V-band  (50-75  GHz),  and  W-band  (75-110  GHz)  fre¬ 
quencies.  The  manufacturer  provided  calibration  curves  in  frequency  and 
power  for  each  head  are  given  in  Figures  16-18.  The  sweeper  system  was 
lent  to  MRC  by  Dan  Hunt  of  White  Sands  Missile  Range  Office  0MEW/P2.  His 
assistance  has  been  invaluable.  Each  of  the  three  millimeter  bands  covered 
by  the  spectrometer  has  its  own  set  of  detectors,  3  each  for  B  and  V-bands, 
and  6  for  W-band.  Type  1N53  crystal  detectors  have  been  chosen  for  all  3 
bands,  due  to  their  low  cost  and  d  availability.  The  crystal  housings 
are  Baytron  #4-22-200  for  B  band,  Baytron  #4-15-2000  for  V  band,  and 
Baytron  #4-10-200  for  W  band.  Each  channel  also  has  its  own  pair  of  grat¬ 
ings  which  have  been  optimized  for  use  at  that  particular  wavelength 
window.  The  spectrometer  operates  in  only  one  band  at  a  time  with  only  one 
set  of  detectors  being  used  at  a  time.  In  each  band  a  nominal  25  dB 
standard  gain  horn  is  used  as  the  input  to  the  spectrometer  in  order  to 
increase  the  signal.  Calibration  curves  for  the  horns  in  each  band  are 
presented  in  Figures  19-21.  Calibration  of  each  band  for  frequency 
response  and  sensitivity  is  carried  out  separately. 
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Figure  16.  Manufacturer  provided  calibration  curve  of  B-band  sweeper  power  and  frequency  output 

as  a  function  of  oscillator  head  anode-cathode  voltage  (Micro-Now  Model  #705B  mainframe 
and  Model  #729  BWO  head).  The  sweeper  drives  the  spectrometer  and  oscilloscope  in  the 
configuration  of  Figures  22  and  23  to  produce  the  calibration  data  shown  in  Figures  24 
through  35. 


Figure  17  Manufacturer  provided  calibration  curve  of  V-band  sweeper  power  and  frequency 
output  as  a  function  of  oscillator  head  anode-cathode  voltage  (Micro-Now  Model 
#7Q5B  mainframe  and  Model  #725  BWO  head). 


Figure  18.  Manufacturer  provided  calibration  curve  of  W-band  sweeper  power  and  frequency 
output  as  a  function  of  oscillator  head  anode-cathode  voltage  (Micro-Now  Model 
#7058  mainframe  and  Model  #728  BWO  head). 
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Manufacturer  provided  calibration  curve  for  the  B-band 
(35-50  GHz)  25  dB  standard  gain  horn  used  in  spectrometer 
calibration  (TRG  Model  #B861 ) . 


Frequency  (GHz) 


Figure  20.  Manufacturer  provided  calibration  curve  for  the  V-band 

(50-75  GHz)  25  dB  standard  gain  horn  used  in  spectrometer 
calibration  (TRG  Model  #V861). 
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Manufacturer  provided  calibration  curve  for  the  W-band 
(75-110  GHz)  25  dB  standard  gain  horn  used  in  spectrometer 
calibration  (TRG  Model  #W861). 
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Two  calibrations  were  performed  for  each  band.  The  flatness  of  the 
frequency  response  of  each  of  the  three  (or  six)  detectors  is  measured 
using  the  apparatus  shown  in  Figure  22.  Two  series  10  or  15  dB  couplers 
are  used  to  attenuate  the  nominal  100  mW  sweeper  output  to  levels  which 
will  not  overdrive  the  crystal.  As  shown  in  the  figure,  the  sweeper  pro¬ 
vides  a  ramp  signal  to  the  oscilloscope  which  varies  linearly  in  time  cor¬ 
responding  to  the  linear  time  variation  of  the  output  frequency.  The  scope 
photo  then  records  the  frequency  response  of  the  crystal,  uncorrected  for 
variations  in  the  sweeper  output  power. 

Next  the  array  of  detectors  is  positioned  inside  the  spectrometer  so 
that  each  samples  a  different  frequency  bin.  In  this  configuration  shown 
in  Figure  23,  the  sweeper  drives  the  spectrometer  directly,  since  the  spec¬ 
trometer  itself  is  very  lossy,  of  order  35  dB.  In  order  to  increase  the 
signal  to  a  measurable  level,  a  25  dB  standard  gain  horn  is  attached  to  the 
waveguide  entrance  flange  which  feeds  each  of  the  crystals.  The  cali¬ 
bration  curves  for  these  horns  are  shown  in  Figures  19,  20,  and  21  for  B, 

V,  and  W-band,  respectively.  Power  variation  in  the  sweeper  output  as  a 
function  of  frequency  is  shown  in  Figures  16,  17,  and  18  for  B,  V,  and 
W-band,  respectively.  To  the  extent  that  the  detector  crystal  response  is 
linear  with  input  power,  these  curves  and  the  following  data  allow  an 
absolute  power  calibration  of  the  spectrometer  frequency  response. 

Raw  data  for  the  three  B-band  detectors'  frequency  response  is  shown 
in  Figure  24.  In  order  to  determine  the  actual  crystal  sensitivity  as  a 
function  of  frequency,  this  raw  data  should  be  divided  point-by-point  by 
the  sweeper  power  output  as  a  function  of  frequency.  The  latter  quantity 
is  available  from  the  plots  of  Figure  16.  Note  that  this  process  also 
involves  an  assumption  that  the  crystal  voltage  output  is  linear  in  the 
power  input,  i.e.,  the  crystal  is  operating  in  the  square  law  regime. 

Since  the  extensive  tests  required  to  verify  square  law  operation  over  the 
entire  frequency  band  have  not  been  done,  the  tedious  point-by-point  divi¬ 
sion  has  not  been  carried  out.  Noting  that  the  input  power  does  not  vary 
by  more  than  a  factor  of  2,  some  broad  conclusions  can  still  be  drawn  from 
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the  raw  data.  First,  detector  response  is  down  considerably  at  high  fre¬ 
quencies,  dropping  by  the  order  of  one  decade  above  44  GHz.  Second, 
response  is  fairly  flat,  within  a  factor  of  2,  in  the  32  to  44  GHz 
frequency  bin.  The  fine  structure  in  the  raw  data  is  due  to  resonances  in 
the  detector  housing  and  some  inherent  non-linearity  in  the  crystal 
response.  Superior  quality  detectors  can  be  purchased  at  a  premium  price 
for  future  work,  if  it  becomes  necessary  to  provide  flatter  response  over 
the  entire  band. 

Next,  detectors  Bl,  B3  and  B2  were  installed  in  the  spectrometer  at 
stations  15.5,  8.0,  and  3.0  inches,  respectively,  and  the  spectrometer 
driven  with  the  sweeper  in  the  configuration  of  Figure  23.  Data  was  taken 
with  the  grating  angle  at  42  degrees  and  48  degrees,  as  shown  in  Figures  25 
and  26.  The  grating  angle  determines  the  acceptance  frequency  of  each 
detector  through  the  optics  criteria  described  in  References  10  and  11. 

A  summary  of  the  B-band  calibration  data  is  presented  in  Table  III  for 
the  42  degree  and  48  degree  grating  angles  tested.  The  center  frequency 
for  peak  sensitivity  and  the  peak  output  voltage  are  listed  for  each  detec¬ 
tor.  Note  that  the  FWHM  response  is  approximately  2  GHz,  yielding  a 
Af/f0  parameter  of  order  0.06.  A  figure  of  merit  parameter  has  been  cal¬ 
culated  by  dividing  the  output  voltage  by  the  sweeper  output  power,  taking 
into  account  the  frequency  dependence  of  the  horn  gain  and  the  loss  due  to 
the  10  dB  coupler  preceding  the  spectrometer.  To  the  extent  that  the 
crystal  output  is  linear  with  input  power,  the  figure  of  merit  is  an 
absolute  calibration  of  the  spectrometer  sensitivity  for  that  particular 
detector,  frequency,  and  grating  angle.  Table  III  shows  that  the  figure  of 
merit  varies  by  nearly  two  decades  from  one  frequency  bin  to  another.  This 
is  just  a  reflection  of  the  poor  flatness  in  the  frequency  response  of  the 
individual  crystals  and  can  be  rectified  by  purchasing  better  housings  and 
better  crystals. 
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Detector  *B1 ,  Angle  42 


32  GHz  50  GHz 


Detector  #B3,  Angle  42: 


32  GHz  50  GHz 


Detector  *B2,  Angle  42° 


32  GHz  50  GHz 


Spectrometer  calibration  data  for  the  three  B-band  detectors 
(detector  B1  at  15.5  inch,  B3  at  8.0  inch,  and  B2  at  3.0  inch) 
taken  at  a  42  degree  grating  angle,  using  the  configuration 
of  Figure  23. 


Detector  *B1,  Angle  48 
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Detector  #B3,  Angle  48° 
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Detector  ?B2,  Angle  48 


32  GHz  50  GHz 


Figure  26.  Spectrometer  calibration  data  for  the  three  B-band  detectors 
(detector  B1  at  15.5  inch,  B3  at  8.0  inch,  and  B2  at  3.0  inch) 
taken  at  a  48  degree  grating  angle. 
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TABLE  III.  SUMMARY  OF  SPECTROMETER  B-BAND  (32-50  GHz)  CALIBRATION  DATA 


PARAMETER 

Detector 

#B1 

Detector 

#B3 

Detector 

#B2 

Detector  Position  (inches) 

15.5 

8.0 

3.0 

For  42°  Grating  Angle: 

Center  Frequency  for  Peak  Sensitivity 
(GHz) 

Output  Voltage  (mV) 

Horn  Gain  (db) 

Input  Power  (mW) 

Calculated  Figure  of  Merit  (Volts/Watt) 

38.4 

180 

24.2 

60 

1.27xl0-2 

43.8 

59 

24.7 

96 

2.31xl0-3 

48.8 

6.2 

25.0 

145 

1.50X10-14 

For  48°  Grating  Angle: 

Center  Frequency  for  Peak  Sensitivity 
(GHz) 

Output  Voltage  (mV) 

Horn  Gain  (db) 

Input  Power  (mW) 

Calculated  Figure  of  Merit  (Volts/Watt) 

33.2 

89 

23.4 

43 

1.05xl0-2 

38.0 

89 

24.1 

57 

6.75xl0-3 

41.4 

134 

24.5 

70 

7.55xl0-3 

Calibration  of  the  V-band  response  of  the  detectors  and  spectrometer 
was  carried  out  analagously  to  the  B-band  calibration  data.  Figure  27 
shows  the  response  of  the  three  individual  crystals  and  Figures  28,  29,  and 
30  show  the  spectrometer  response  for  44,  46,  and  49  degree  grating  angles, 
respectively.  Peak  frequency  response  and  figure  of  merit  numbers  are  sum¬ 
marized  in  Table  IV.  Note  that  the  separate  detector  responses  are  far 
from  flat,  typically  varying  by  the  order  of  one  decade  between  peaks  and 
valleys.  This  is  the  main  culprit  in  the  factor  of  40  variation  in  the 
sensitivity  of  the  various  frequency  bins.  The  data  shows  typical  FWHM  of 
the  frequency  response  is  about  2  GHz,  yielding  a  Af/f0  parameter  of 
approximately  0.04. 

Calibration  of  the  W-band  detectors  was  carried  out  in  a  similar  man¬ 
ner  to  the  B  and  V-bands,  but  was  complicated  by  the  extremely  hashy  and 
non-uniform  frequency  response  of  these  crystals  and  their  housings.  The 
data  of  Figures  31A  and  318  shows  a  variation  of  more  than  two  decades  in 
the  raw  frequency  response  data  over  the  75-110  GHz  frequency  band  for 
several  of  the  6  detectors.  Calibration  data  for  the  W-band  spectrometer 
configuration  is  shown  in  Figures  32  through  35,  with  a  data  summary  in 
Table  V.  Frequency  response  FWHM  is  typically  2  GHz,  yielding  a  Af/f0 
parameter  of  order  0.027.  The  figure  of  merit  parameter  continues  to  show 
large  variation,  roughly  two  decades,  over  the  frequency  band.  Thus,  the 
frequency  selectivity  of  this  spectrometer  is  good  although  it  displays 
wide  variation  in  power  sensitivity.  The  latter  undesirable  characteristic 
can  be  corrected  in  future  work  by  installing  crystals  with  flatter  fre¬ 
quency  response  in  housings  with  weaker  resonance  characteristics.  How¬ 
ever,  in  its  present  form  and  with  the  calibration  data  presented  here,  the 
spectrometer  is  a  unique  diagnostic  for  characterizing  the  time  dependence, 
power  and  spectral  content  of  pulsed,  millimeter  wavelength  microwave 
sources. 
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Detector  *V2,  Angle  44° 
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Detector  *V3,  Angle  44 
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Figure  28.  Spectrometer  calibration  data  for  the  three  V-band  detectors 
(detector  VI  at  15.5  inch,  V2  at  8.0  inch,  and  V 3  at  3.0  inch) 
taken  at  a  44  degree  grating  angle,  using  the  configuration 
of  Figure  23. 
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Figure  29.  Spectrometer  calibration  data  for  the  three  V-band  detectors 
(detector  VI  at  15.5  inch,  V2  at  8.0  inch,  and  V 3  at  3.0  inch) 
taken  at  a  46  degree  grating  angle. 
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Detector  *V1,  Angle  49° 


50  GHz  75  GHz 


Detector  tfVc,  Angle  49' 


50  GHz  75  GHz 


Detector  *V 3,  Angle  49J 


Figure  30.  Spectrometer  calibration  data  for  the  three  V-band  detectors 
(detector  VI  at  15.5  inch,  V2  at  8.0  inch,  and  V 3  at  3.0  inch) 
taken  at  a  49  degree  grating  angle. 
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TABLE  IV.  SUMMARY  OF  SPECTROMETER  V-BAND  (50-75  GHz)  CALIBRATION  DATA 


Detector 

#V1 

Detector 

#V2 

Detector 

#V3 

Detector  Position  (inches) 

15.5 

8.0 

3.0 

For  44°  Grating  Angle: 

Center  Frequency  for  Peak 
Sensitivity  (GHz) 

59.5 

67.9 

73.0 

Output  Voltage  (mV) 

20 

3.5 

1.2 

Horn  Gain  (db) 

24.8 

25.1 

25.2 

Input  Power  (mW) 

82 

127 

no 

Figure  of  Merit  (Volts/Watt) 

8.97xl0-4 

9.46xl0-5 

3.05xl0-5 

For  46°  Grating  Angle 

Center  Frequency  for 

Peak  Sensitivity  (GHz) 

57.1 

63.8 

70.2 

Output  Voltage  (mV) 

11.5 

3.8 

4.3 

Horn  Gain  (db) 

24.5 

24.9 

25.1 

Input  Power  (mW) 

66 

96 

126 

Figure  of  Merit  (Volts/W) 

6.87xl0-4 

1.42xl0-4 

1.17xl0-4 

For  49°  Grating  Angle 

Center  Frequency  for 

Peak  Sensitivity  (GHz) 

53.1 

60.2 

64.7 

Output  Voltage  (mV) 

15.8 

13.8 

9.6 

Horn  Gain  (db) 

24.2 

24.7 

25.0 

Input  Power  (mW) 

50 

80 

105 

Figure  of  Merit  (Volts/Watt) 

1.33xlO'3 

6.49xl0-4 

3.21xl0-4 

t 

75  GHz 
Figure  31A. 
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Detector  calibration  data  for  the  first  three  W-band  detectors 
W1 ,  W2,  and  W3,  taken  using  the  apparatus  of  Figure  22.  Data 
shown  on  the  right  is  taken  at  several  times  the  gain  of  the  left, 
as  is  reflected  in  the  labels  on  the  scope  photos. 


60 


75  GHz 
Figure  31B. 


110  GHz 


75  GHz 


110  GHz 


Detector  calibration  data  for  the  last  three  W-band  detectors,  W4 , 
W5.  and  W6,  taken  using  the  apparatus  of  Figure  22.  Data  shown 
on  t!  a  right  is  taken  at  several  times  the  gain  of  the  left,  as  is 
reflected  in  the  labels  on  the  scope  photos. 


Detector  #W2,  Angle  44: 
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Detector  *W3,  Angle  44 


Detector  *W6,  Noise  Only 
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110  GHz 


75  GHz 


110  GHz 


Figure  32.  Spectrometer  calibration  data  for  the  W-band  detectors  W2,  W3,  and 
W6  (detector  W2  at  5.5  inch,  W3  at  8.0  inch,  and  W6  at  15.5  inch) 
taken  at  a  44  degree  grating  angle,  using  the  configuration  of  Figure  23 
Left  photo  is  signal  plus  noise,  right  photo  is  noise  only,  taken  with 
microwave  absorbing  material  covering  the  detector. 
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Figure  33B.  Spectrometer  calibration  data  for  the  W-band  detectors  W1-W6  (detectors 
W1-W6  located  at  3.0,  5.5,  8.0,  10.2,  13.0,  and  15.5  inch,  respectively) 
taken  at  a  46  degree  grating  angle.  Left  photo  is  signal  plus  noise, 
right  photo  is  noise  only,  taken  with  microwave  absorbing  material 
covering  the  detector. 
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Figure  34A.  Spectrometer  calibration  data  for  the  W-band  detectors  W1-W6  (detectors 
W1-W6  located  at  3.0,  5.5,  8.0,  10.2,  13.0,  and  15.5  inch,  respectively 
taken  at  a  49  degree  grating  angle.  Left  photo  is  signal  plus  noise, 
right  photo  is  noise  only,  taken  with  microwave  absorbing  material 
covering  the  detector. 
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Figure  34B.  Spectrometer  calibration  data  for  the  W-band  detectors  W1-W6  (detectors 
W1-W6  located  at  3.0,  5.5,  8.0,  10.2,  13.0,  and  15.5  inch,  respectively) 
taken  at  a  49  degree  grating  angle.  Left  photo  is  signal  plus  noise, 
right  photo  is  noise  only,  taken  with  microwave  absorbing  material 
covering  the  detector. 
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Figure  35A.  Spectrometer  calibration  data  for  the  W-band  detectors  W1-W6  {detectors 
W1-W6  located  at  3.0,  5.5,  8.0,  10.2,  13.0,  and  15.5  inch,  respectively) 
taken  at  a  51.5  degree  grating  angle.  Left  photo  is  signal  plus  noise, 
right  photo  is  noise  only,  taken  with  microwave  absorbing  material 
covering  the  detector. 
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35B .  Spectrometer  calibration  data  for  the  W-band  detectors  W1-W6  (detectors 
W1-W6  located  at  3.0,  5.5,  8.0,  10.2,  13.0,  and  15.5  inch,  respectively) 
taken  at  a  51.5  degree  grating  angle.  Left  photo  is  signal  plus  noise, 
right  photo  is  noise  only,  taken  with  microwave  absorbing  material 
covering  the  detector. 


TABLE  V.  SUMMARY  OF  SPECTROMETER  W-BAND  (75-110  GHz)  CALIBRATION  DATA 


Calculated  Figure  of  Merit  (Volts/Watt)  2.7 6x10”  5.52xlO“4  2.12xl0"4  5.58xl0”5  I  1.87xlO“5  5.22x10” 


SECTION  IV 

ELECTRON  BEAM  SOURCE 


The  design  of  the  pulse  power  system  which  drives  the  vircator  e-beam 
diode  has  been  previously  described.11  A  schematic  of  the  electrical 
system  is  shown  in  Figure  36.  Briefly,  the  electrical  equivalent  of  a 
Blumlein  type  pulse  forming  line  (PFL)  is  built  from  discrete  capacitors 
and  inductors.  The  Blumlein  is  DC  charged  to  nominally  50  kV,  and  the  PFL 
is  erected  by  a  triggered  spark  gap  switch  which  is  driven  by  a  time  delay 
generator.  The  purpose  of  the  time  delay  generator  is  to  trigger  the 
e-beam  near  peak  magnetic  field  time,  which  is  typically  600  ysec  after  the 
magnet  circuit  is  fired.  The  artificial  Blumlein,  as  built  inside  a  high 
voltage  insulated  screen  box,  together  with  the  charging  and  control  com¬ 
ponents,  is  shown  in  Photo  4.  A  3  cm  diameter,  50  ohm  high  voltage  cable 
connects  the  Blumlein  source  output  to  the  vircator  vacuum  vessel  shown  in 
Photo  5. 


An  important  consideration  in  the  vircator  diode  development  work  has 
been  to  find  a  material  which  exhibits  fast  electron  field  emission  turn  on 
and  slow  gap  closure.  The  gap  closure  problem  is  severe  because  the  low 
vircator  voltage  requires  operation  at  very  small  anode-cathode  gaps,  of 
order  a  few  mm,  in  order  to  reach  the  threshold  electric  fields  for  field 
emission.  Carbon  has  been  tried  as  the  cathode  material  because  of  its 
well-known  low  threshold  for  field  emission,  of  order  50  kV/cm.  However, 
the  resulting  slow  turn  on,  spotty  emission,  and  high  gap  closure  rate 
(several  cm/ysec)  made  carbon  unacceptable  for  the  cathode  material. 

An  acceptable  cathode  material  developed  in  the  course  of  this  work  is 
the  cloth  fiber  array.  It  has  shown  fast  turn-on,  as  well  as,  very  slow 
gap  closure  rates,  less  than  1  cm/ysec.  A  typical  voltage  and  current 
trace  for  the  cloth  fiber  array  is  shown  in  Figure  37.  This  data  was  taken 


ELECTRON  BEAM  PULSED  POWER  SUPPLY 


0 


Figure  36.  Schematic  of  the  limped  element  Blumlein  pulse- formi ng  line.  Eight  stages 
are  used  to  give  a  pulse  length  of  600  nsec  at  a  50  ohm  impedance. 
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Photograph  4. 


The  8  stage  artificial  Blumlein  is  shown  inside  the  high 
voltage  insulated  screen  box,  together  with  charging  and 
triggering  circuitry.  The  16  inductors,  oriented  vertically 
are  clearly  visible.  The  48  2700  pf,  40  kV  rated  button 
capacitors,  stacked  in  sets  of  3  each  located  under  the 
white  polyethylene  insulation,  are  not  visible. 
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the  e-beam  pulser  to  the  vacuum 
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Figure  37.  Voltage  and  current  traces  for  a  50  kV  charge  voltage  and  a  4  mm  A-K  gap. 

Voltage  risetime  is  20  nsec.  Peak  voltage  and  currents  are  45  kV  and  1.07  kA, 
respectively. 
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with  a  1.9  cm  diameter-cathode,  a  grid  type  anode,  and  a  4  mm  anode-cathode 
gap  operating  as  a  standard  planar  diode.  Emission  is  within  a  factor  of 
20%  of  the  space-charge  limit.  Note  that  the  current  turns  on  as  fast  as 
the  voltage,  approximately  20  nanoseconds,  and  that  the  gap  stays  open  for 
essentially  the  entire  pulse  length,  0.6  usee.  These  characteristics  make 
the  cloth  fiber  array  the  favored  choice  for  the  vircator  cathode  relative 
to  the  other  cold  cathode  emitters. 

Examination  of  the  cloth  fiber  array  with  an  optical  microscope  shows 
that  it  is  composed  of  a  regular  array  of  approximately  3600  nonconducting 
thin  fibers  per  cm2  with  the  fibers  oriented  perpendicular  to  a  partially 
transparent  nonconducting  cloth  backing.12  The  emission  region  then  con¬ 
sists  of  an  array  of  fibers  perpendicular  to  a  conducting  cathode  surface. 

A  surface  flashover  along  the  individual  fibers  contributes  to  the  plasma 
formation  process,  and  this  apparently  results  in  the  formation  of  more 
emission  sites  than  for  the  analogous  carbon  cold  cathode.  Both  long  life¬ 
time  and  the  low  gap  closure  velocity  are  the  result  of  reduced  energy  per 
emission  site.  The  amount  of  energy  per  site  is  very  small,  because  the 
breakdown  mechanism  is  effectively  a  capacitive  discharge  of  low  capaci¬ 
tance.  In  this  sense  we  believe  that  the  cloth  fiber  array  acts  like  the 
Corona  electron  source  developed  by  Hel ionetics13  for  laser  pre-ioniza¬ 
tion.  The  axial  surface  flashover  mechanism  is  more  desirable  than  the 
conventional  cold  cathode  emission  process,  because  production  of  plasma  in 
this  manner  inhibits  the  formation  of  surface  cathode  spots. 
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APPENDIX  A 


LISTING  OF  FORTRAN  SOURCE  CODE  FOR 
MAGNETIC  FIELD  CALCULATION 


SOLENOIDAL  MAGNETIC  FIELD  CALCULATION 

BASE  CODE  BY  JOHN  FREEMAN,  MODIFIED  BY  DON  VOSS  FOR  FIELD  LINE 
CALCULATION. 

COMMON  /A/  RB( 50) ,ZB(50) ,RF(50) ,ZF ( 50) 

COMMON  /B/  RCC (100) ,ZCL ( 100)  ,ZCR( 100) , X >  *  T( 100) ,CCC( 100) 

COMMON  / V0SS1/  BRAD(41,  201),  BAX IAL(41 ,  201) 

COMMON  / VOS S2 /  RADIUS,  XLENG,  DR,  DZ,  DRXDZ 
DIMENSION  I  HEAD ( 7 ) 

DIMENSION  STO ( 5000 ) 

DIMENSION  XPLOT ( 1001 ) ,  YPLOT(lOOl) 

COMMON  /C/  R ( 2001  ,Z ( 200) ,NR ,NZ ,NCT 
DIMENSION  IPHEAD(3) 

REAL  LENGTH 

DATA  IPHEAD  /4HBRJ- .4HBZJ- ,4HPSJ=/ 

OPEN  (UNIT-11,  NAME- 1 MAGIN.ECH1 ,  STATUS-' NEW' ) 

OPEN  (UNIT-12,  NAME-'MAGOUT.LPT' ,  STATUS-' NEW ) 

PI  =  3.1415926535 
RADDEG  =  180.  /  PI 

PRINT  * , 1  ALL  DIMENSIONS  ARE  IN  CENTIMETERS' 

100  FORMAT ( 7A10) 

PRINT  ENTER  HEADING  FOR  THIS  PROBLEM' 

READ100.IHEAD 

101  FORMAT ( IX ,7A10/ / ) 

PRINT  *, ' FOR  KEYBOARD  INPUT  TYPE  1;  OTHERWISE  TYPE  O' 

READ  * , I NP 

PRINT  *, ' FOR  SHORT  PRINT  TYPE  1;0THERWISE  TYPE  0' 

READ  *,KPT 

IF  ( INP.EQ.O)  GO  TO  30 

PRINT  *,' ENTER  RADIUS  OF  COMPUTATION’ 

READ  *, RADIUS 

PRINT  *, 'ENTER  LENGTH  OF  COMPUTATION' 

READ  *, LENGTH 

PRINT  ENTER  NUMBER  OF  RADIAL  MESH  LINES' 

READ  *,NR 

PRINT  *,'  ENTER  NUMBER  OF  AXIAL  MESHLINES' 

READ  *,NZ 
NCT-0 

2  PRINT  *, 'ENTER  COIL  IDENT  NUMBER;  0  TO  TERMINATE' 

READ  *,N C 


f 


PRINT  ENTER  NUMBER  OF  TURNS' 

READ  *,  XNTC(NC) 

PRINT  *, 'ENTER  COIL  CURRENT  IN  AMPERES' 

READ  *,CCC(NC) 

GO  TO  2 

30  CONTINUE 

OPEN  (UNIT=1 ,  NAME= 'MAGIN.DAT' ,  STATUS= 'OLD* ) 

READ(1,*)  RADIUS 
READ(1,*)  LENGTH 
READ ( 1 ,* )  NR 
READ( 1,*)  NZ 
READ ( 1 , * )  NCT 
%  DO  32  1=1,  NCT 

READ( 1,*)  RCC(I),  ZCL(I),  ZCR(I),  XNTC(I),  CCC(I) 

32  CONTINUE 
1  CONTINUE 

400  FORMAT  (IX,  1PE16.8) 

401  FORMAT  (IX,  16) 

402  FORMAT  (IX,  5( 1PE16.8) ) 

WRITE  (11,400)  RADIUS 
WRITE  (11,400)  LENGTH 
WRITE  (11,401)  NR 
WRITE  (11,401)  NZ 
WRITE  (11,401)  NCT 

#  DO  31  1=1,  NCT 

WRITE  (11,402)  (RCC(I),  ZCL(I),  ZCR(I),  XNTC(I),  CCC(I)) 

31  CONTINUE 
XLENG  =  LENGTH 
DR  =  0. 

DZ  =  0, 

#  IF  (NR  .EQ.  1)  GO  TO  9 
NR1  =  NR  -  1 

DR  =  RADIUS  /  NR1 
R( 1)  =  0. 

DO  3  1=2,  NR 

3  R ( I )  =  R(I-l)  +  DR 

,  9  IF  (NZ  .EQ.  1)  GO  TO  8 

NZ1  =  NZ  -  1 
DZ  =  LENGTH  /  NZ1 
Z  ( 1 )  =  0. 

DO  4  1=2,  NZ 

4  Z(I)  =  Z(I-l)  +  DZ 

#  8  CONTINUE 

DRXDZ  =  DR  *  DZ 
PRINT  201 

201  FORMAT  ('  COMPUTATION  REGION  DATA') 

PRINT  102, RADIUS, LENGTH, DR, DZ, NR, NZ 
102  FORMAT  (/,'  RADIUS=',  1PE12.3,  5X,  'LENGTH= ' ,  1PE12.3,  4X, 


9 


202  FORMAT  ('  COIL  DATA') 

00  5  1=1,  NCT 

PRINT  103,  (RCC(I),  ZCL(I),  ZCR(I),  XNTC(I),  CCC(I)) 

103  FORMAT  ('  RCC=',  1PE9.2,  '  ZCL=',  1PE9.2, 

1  '  ZCR= 1 ,  1PE9.2,  '  NTC= 1 ,  1PE9.2,  '  CCC=',  1PE9.2) 

5  CONTINUE 

IF  (KPT  .EQ.  1)  GO  TO  33 
PRINT  105 

105  FORMAT (/,  '  RADIAL  MESH  LINE  LOCATIONS') 

PRINT  104,  (R( I ),  1=1,  NR) 

PRINT  106 

106  FORMAT  (/,  '  AXIAL  MESH  LINE  LOCATIONS') 

PRINT  104,  (Z(I),  1=1,  NZ) 

104  FORMAT  ( ( 1 0 ( 5 ( 1 PE  1 3 . 5 ) / ) )/ ) 

33  CONTINUE 

WRITE  (12,  201) 

WRITE  (12,  102)  RADIUS,  LENGTH,  DR,  DZ,  NR,  NZ 
WRITE  (12,  202) 

DO  250  1=1,  NCT 

250  WRITE  (12,  103)  RCC(I),  ZCL(I),  ZCR(I),  XNTC(I),  CCC(I) 

WRITE  (12,  105) 

WRITE  (12,  104)  (R(  I ),  1  =  1,  NR) 

WRITE  (12,  106) 

WRITE  (12,  104)  (Z(I),  1=1,  NZ) 

NT  =  NR  *  NZ 

PRINT  *,  '  ENTER  -2  TO  BYPASS,  -1  FOR  BR,  0  FOR  BZ,  1  FOR  PS  I ' 
READ  *,  I F I  ELD 

IF  (IFIELD  .LE.  -2)  GO  TO  16 
CALL  FIELD  (STO,  IFIElO) 

17  PRINT  *,  '  ENTER  J-LINE  FOR  PRINT;  0  TO  TERMINATE' 

IPT  =  IFIELD  +  2 
READ  *,  J 

IF  (J  .EQ.  0)  GO  TO  16 
PRINT  108,  IPHEAD(IPT),  J 
108  FORMAT  (//IX,  A7,  13) 

DO  7  L=1 ,  NZ 
JL  =  (L-l)  *  NR  +  J 
7  RB(L)  =  STO(JL) 

PRINT  104,  (RB(L) ,  L=l,  NZ) 

GO  TO  17 
16  CONTINUE 

WRITE  (6,  300) 

300  FORMAT  (/,  '  TYPE  A  1  IF  YOU  WANT  TO  DO  FIELD  LINES') 

READ  (5,  301)  IFLINE 

301  FORMAT  (12) 

IF  (IFLINE  .EQ.  0)  GO  TO  399 
IFIELD  =  -1 

CALL  FIELD  (STO,  IFIELD) 


JL  =  (L  -  1)  *  NR  +  K 
BRAD(K,L)  =  STO(JL) 

364  CONTINUE 

365  CONTINUE 
IF  I  ELD  =  0 

CALL  FIELD  (STO,  I F I  ELD ) 

DO  367  K=l,  NR 
DO  366  L=l,  NZ 
JL  =  (L  -  1)  *  NR  +  K 
BAXIAL(K.L)  =  STO(JL) 

366  CONTINUE 

367  CONTINUE 
WRITE  (12,  305) 

305  FORMAT  (//,  '  LISTING  OF  BR  AND  BZ  OVER  THE  COMPUTATIONAL  GRID') 
DO  369  K=l,  NR 

WRITE  (12,  306) 

306  FORMAT  (/,'  K\  3X,  * L * ,  5X,  'R',  9X,  *Z * ,  6X,  *BR(K,L)  1 ,  4X, 

1  ’BZ(K,L)  ' ,  4X,  'B(K,L) ' ,  5X,  'BR/BZ',  5X,  'THETA') 

XKV  =  K  -  1 
RADV  =  XKV  *  DR 
DO  368  L=1 ,  NZ 
XLV  =  L  -  1 
ZV  =  XLV  *  DZ 
RCOMP  =  BRAD(K,L ) 

ZCOMP  =  BAXIAL(K.L) 

BMAG  =  SQRT  (RCOMP  *  RCOMP  +  ZCOMP  *  ZCOMP) 

BRATIO  =  0. 

IF  (ZCOMP  .NE.  0.)  BRATIO  =  RCOMP  /  ZCOMP 
ANGV  =  RADDEG  *  ATAN(BRATIO) 

WRITE  (12,  307)  K,  L,  RADV,  ZV,  RCOMP,  ZCOMP,  BMAG,  BRATIO,  ANGV 

307  FORMAT  (IX,  12,  14,  F9.5,  F9.4,  3F11.3,  IX,  1PE10.3,  0PF9.4) 

368  CONTINUE 

369  CONTINUE 

READ  DATA  FOR  FIELD  LINE  CALCULATION 
380  WRITE  (6,  500) 

500  FORMAT  ('  ENTER  A  1  FOR  FIELD  LINE  TTY  INPUT,  0  FOR  DISK  INPUT,', 

1  '  -1  TO  TERM  —  ',  $) 

READ  (5,  *)  IDVOPT 
IF  (IDVOPT  .EQ.  -1)  GO  TO  590 
IF  (IDVOPT  .EQ.  0)  GO  TO  555 
550  WRITE  (6,  501) 

501  FORMAT  ('  HOW  MANY  FIELD  LINES  (MAX  =  20)  --  ',  $) 

READ  (5,  *)  IFLNUM 

IF  (IFLNUM  .LE.  0)  GO  TO  550 
552  WRITE  (6,  502) 

502  FORMAT  ('  ENTER  FIELD  LINE  STEP  DISTANCE  DS  (CM)  —  ',  $) 


553  WRITE  (6,  503) 

503  FORMAT  ('  ENTER  PLOT  SPACING  DISTANCE  DPLOT  (CM)  —  ',  $) 

READ  (5,  *)  DPLOT 

IF  (DPLOT  .LT.  0.)  GO  TO  553 

554  WRITE  (6,  504) 

504  FORMAT  ('  ENTER  RADIAL  RANGE  (2  NUMBERS)  FOR  FIELD  LINES  (CM)  — 

1  ,  $) 

READ  (5,  *)  RFLBEG,  RFLEND 
IF  (RFLBEG  .GT.  RFLEND)  GO  TO  554 
567  WRITE  (6,  517) 

517  FORMAT  ('  ENTER  RADIAL  GRAPHICS  ENHANCEMENT  FACTOR  —  $) 

READ  (5,  *)  RFACT 

IF  (RFACT  .LE.  0.)  GO  TO  567 

WRITE  (6,  505) 

505  FORMAT  ('  ENTER  R  AND  Z  DISTANCES  FOR  PLOT  WINDOW', 

1  '  (2  NUMBERS)  —  $) 

READ  (5,  *)  RPLOTW,  ZPLOTW 
GO  TO  556 

555  CONTINUE 

READ  (1,  *)  IFLNUM 

READ  (1,  *)  DS 

READ  (1,  *)  DPLOT 

READ  (1,  *)  RFLBEG,  RFLEND 

READ  (1,  *)  RFACT 

READ  (1,  *)  RPLOTW,  ZPLOTW 

556  CONTINUE 

WRITE  (12,  510)  IFLNUM,  DS,  RFLBEG,  RFLEND,  DPLOT 
510  FORMAT  (//,  '  FIELD  LINE  PLOT  SPECIFICATIONS  AS  FOLLOWS:',  /, 

1  IX,  12,  '  FIELD  LINES  CALCULATED  WITH  DS  =',  1PE14.7, 

2  '  CM,',/,  '  BEGINNING  AT  A  RADIUS  OF*,  0PF10.6,  '  CM  AND 

3  /,  '  ENDING  AT  A  RADIUS  OF',  0PF10.6,  '  CM.  FIELD  LINE' 

4  /,  '  (R,Z)  POSITIONS  RECORDED  AT  INTERVALS  OF',  0PF10.6, 

5  'CM.') 

IF  (RPLOTW  .GT.  0.)  GO  TO  570 
RPLOTW  =  RADIUS 
DO  569  IDUM=1 ,  NCT 

IF  (RCC(IDUM)  .GT.  RPLOTW)  RPLOTW  =  RCC(IDUM) 

569  CONTINUE 

570  IF  (ZPLOTW  .GT.  0.)  GO  TO  572 
ZPLOTW  =  XLENG 

DO  571  IDUM=1 ,  NCT 

IF  (ZCR(IDUM)  .GT.  ZPLOTW)  ZPLOTW  =  ZCR(IDUM) 

571  CONTINUE 

572  WRITE  (12,  523)  RFACT,  RPLOTW,  ZPLOTW 

523  FORMAT  (’  RADIAL  SCALE  ON  PLOTS  ENHANCED  BY',  1PE14.7,  /, 

1  '  (R, Z)  PLOT  RANGE  =  (',  0PF9.5,  0PF9.5,  ').') 

WRITE  (6,  510)  IFLNUM,  DS,  RFLBEG,  RFLEND,  DPLOT 
WRITE  (6,  523)  RFACT,  RPLOTW,  ZPLOTW 


WRITE  (11,  400)  DPLOT 

WRITE  (11,  402)  RFLBEG,  RFLEND 

WRITE  (11,  400)  REACT 

WRITE  (11,  402)  RPLOTW,  ZPLOTW 

START  FIELD  LINES  AT  Z=0  AND  USE  'SHOOTING*  TECHNIQUE  FOR  INTEGRATING 


DRFL  =  0. 

XFLNUM  =  IFLNUM 

IF  (IFLNUM  .GT.  1)  DRFL  =  (RFLEND  -  RFLBEG)  /  (XFLNUM  -  1.) 

PLOT  SET-UP  —  PLOT  AXES  AND  LABELS  —  11-30-82 

XPLOT(l)  =  0. 

YPLOT(l)  =  0. 

XPL0T(2)  =  ZPLOTW 
YPLOT (2)  =  RPLOTW 
RADEFF  =  RPLOTW  *  RFACT 
CALL  STARTT ( 1) 

CALL  PLOJB  (XPLOT,  YPLOT,  2,  1,  0,  0,  0.0,  ZPLOTW, 

1  RADEFF, 37H  VIRCATOR  VACUUM  FIELD  LINE  STRUCTURE, 

2  -37.25H  AXIAL  CO-ORDINATE  Z  (CM), 

3  25,  2H  R,  2) 

DO  560  IDUM=1 ,  IFLNUM 

XDUM  *  IDUM  -  1 

RSTART  =  RFLBEG  +  XDUM  *  DRFL 

ZSTART  =  0. 

WRITE  (12,  511)  IDUM 

511  FORMAT  (//,  *  FIELD  LINE  NUMBER',  14,  '  (R,Z)  TRAJECTORY  IS  AS' 

1  '  FOLLOWS  :') 

WRITE  (12,  512) 

512  FORMAT  (/,'  STEP  #',  5X,  'R\  9X,  'Z',  6X,  'BR(K.L) ' ,  4X, 

1  'BZ(K.L) ' ,  4X,  'B(K,L) ',  5X,  'BR/BZ* ,  5X,  'THETA') 

DISTDV  =  DPLOT 
INC  =  0 

DO  559  IDV=1 ,  10001 

CALL  BINTRP  (RSTART,  ZSTART,  BRINT,  BZINT,  IOFF) 

IF  (IOFF  .EQ.  1)  GO  TO  559 

BMAG  =  SQRT  (BRINT  *  BRINT  +  BZINT  *  BZINT) 

IF  (DISTDV  .LT.  DPLOT)  GO  TO  566 
BRATIO  =  0. 

IF  (BZINT  .NE.  0.)  BRATIO  =  BRINT  /  BZINT 
ANGV  =  RADDEG  *  ATAN  (BRATIO) 

WRITE  (12,  513)  IDV,  RSTART,  ZSTART,  BRINT,  BZINT,  BMAG,  BRATIO, 
1  ANGV 

513  FORMAT  (IX,  16,  F9.5,  F9.4,  3F11.3,  IX,  1PE10.3,  0PF9.4) 

DISTDV  =  DISTDV  -  DPLOT 


YPLOT (INC)  =  RSTART 
566  IF  (BMAG  .GT.  0.)  GO  TO  557 
WRITE  (12,  515) 

515  FORMAT  ('  RUN  TERMINATED  DUE  TO  ZERO  B  FIELD') 

GO  TO  561 

557  FACDV  =  DS  /  BMAG 

RSTART  =  RSTART  +  BRINT  *  FACDV 
IF  (RSTART  .GT.  0.)  GO  TO  558 
WRITE  (12,  514) 

514  FORMAT  ('  RUN  TERMINATED  DUE  TO  PROJECTION  THROUGH  ZERO  RADIUS') 
GO  TO  561 

558  ZSTART  =  ZSTART  +  BZINT  *  FACDV 
DISTDV  =  DISTDV  +  DS 

559  CONTINUE 

561  IF  (INC  .LE.  1)  GO  TO  560 

PLOT  THE  FIELD  LINE  JUST  CALCULATED 

CALL  CHOPAR  (XPLOT,  YPLOT,  INC,  ZPLOTW,  RPLOTW) 

CALL  PLOJB  (XPLOT,  YPLOT,  INC,  1,  0,  -1,  0.0,  ZPLOTW,  RADEFF , 

1  0,  0,  0,  0,  0,  0) 

560  CONTINUE 
399  CONTINUE 

PLOT  THE  POSITION  OF  EACH  TURN  OF  THE  COIL  SET 

DO  575  IDV=1,  NCT 
INC  =  XNTC(IDV) 

XPLOT ( 1)  =  ZCL(IDV) 

YPLOT(l)  =  RCC(IDV) 

IF  (INC  .EQ.  1)  GO  TO  574 

DXCOIL  =  (ZCR(IDV)  -  ZCL(IDV))  /  (XNTC(IDV)  -  1.) 

DO  573  IDUM=2,  INC 
XDUM1  =  IDUM  -  1 

XPLOT ( IDUM)  =  ZCL(IDV)  +  XDUM1  *  DXCOIL 

573  YPLOT ( IDUM)  =  RCC(IDV) 

574  ICHARP  =  -  IDV 

CALL  CHOPAR  (XPLOT,  YPLOT,  INC,  ZPLOTW,  RPLOTW) 

CALL  PLOJB  (XPLOT,  YPLOT,  INC,  1,  -1,  ICHARP,  0.0, 

1  ZPLOTW,  RADEFF,  0,  0,  0,  0,  0,  0) 

575  CONTINUE 

MAKE  HARDCOPY  OF  PLOT  AND  RELEASE  PLOTTER 

CALL  HDCOPY 
CALL  STOPT 


READ  *,  ICORR 

IF  (ICORR  .EQ.  0)  GO  TO  20 

PRINT  *, 'FIELD  CORRECTION  FOR  BZ  ONLY' 

CALL  ECORR(STO) 

CALL  FIELD(STO,IFIELD) 

21  PRINT  VENTER  J-LINE  FOR  PRINT;  0  TO  TERMINATE' 
READ  *,J 

IF(J.EQ.O)  GO  TO  22 
PRINT  108,IPHEAD(IPT),J 
DO  70  L=1 ,NZ 
JL=(L-1 )*NR+J 
70  RR (L )=STO( JL ) 

PRINT  104, (RB(L )  ,L  =  1  ,NZ ) 

GO  TO  21 

22  CONTINUE 
20  CONTINUE 

STOP 

END 


SUBROUTINE  TO  CHOP  ARRAYS 

SUBROUTINE  CHOPAR  (XPLOT,  YPLOT,  NUM,  XMAX,  YMAX ) 
DIMENSION  XPLOT ( 1001 ) ,  YPLOT(lOOl) 

DO  1  1=1,  NUM 

IF  (XPLOT ( I )  .GT.  XMAX)  XPLOT(I)  =  XMAX 
IF  (YPLOT ( I )  .GT.  YMAX)  YPLOT(I)  =  YMAX 
1  CONTINUE 
RETURN 
END 


SUBROUTINE  TO  DO  2-DIMENSIONAL  MAGNETIC  FIELD  INTERPOLATION 
D.  VOSS  -  11/24/82. 

SUBROUTINE  BINTRP  (RCOORD,  ZCOORD,  BRINT,  BZ I  NT ,  IOFF) 

COMMON  /V0SS1/  BRAD(41 ,  201),  BAXIAL (41 ,  201) 

COMMON  /V0SS2/  RADIUS,  XLENG,  DR,  DZ,  DRXDZ 
COMMON  /C/  R(200) ,  Z(200),  NR,  NZ,  NCT 
KUPR  =  (RCOORD  /  DR)  +  2. 

LUPR  =  (ZCOORD  /  DZ)  +  2. 

IF  (KUPR  .LE.  NR  .AND.  LUPR  .LE.  NZ)  GO  TO  3 

1  BRINT  =  0. 

BZINT  =  0. 

IOFF  =  1 

IOFF  =  1  INDICATES  THAT  THE  (R,Z)  POINT  REQUESTED  IS  NOT  ON  THE  GRID 

2  RETURN 

3  IF  (KUPR  .LE.  1  .OR.  LUPR  .LE.  1)  GO  TO  1 


CALCULATE  AREAL  WEIGHTS  -  THIS  SCHEME  ACCURATE  TO  2ND  ORDER  IN  DR,  DZ 

RDIST  =  R{KUPR )  -  RCOORD 
ZDIST  =  Z(LUPR)  -  ZCOORD 
DRDIST  *  DR  -  RDIST 
DZDIST  =  DZ  -  ZDIST 
WLL  =  RDIST  *  ZDIST 
WLR  =  RDIST  *  DZDIST 
WUL  =  DRDIST  *  ZDIST 
WUR  =  DRDIST  *  DZDIST 
KUPRl  =  KUPR  -  1 
LUPRl  =  LUPR  -  1 

BRINT  =  (BRAD(KUPR1, LUPRl)  *  WLL  +  BRAD(KUPR1 ,LUPR )  *  WLR  + 

1  BRAD(KUPR, LUPRl)  *  WUL  +  BRAD (KUPR, LUPR)  *  WUR)  /  DRXDZ 

BZINT  =  (BAXIAL(KUPR1, LUPRl)  *  WLL  +  BAXIAL (KUPRl , LUPR )  *  WLR  + 

1  BAXIAL(KUPR, LUPRl)  *  WUL  +  BAXIAL(KUPR.LUPR)  *  WUR)/DRXDZ 

GO  TO  2 
END 


SUBROUTINE  FIELD(STO, IFIELD ) 

COMMON  /B/  RCC(IOO) ,  ZCL(IOO),  ZCR(IOO),  XNTC(IOO),  CCC(IOO) 
COMMON  /C/  R(200),  Z(200),  NR,  NZ,  NCT 
DIMENSION  STO(l) 

REAL  K,K2 
TPI=6. 2831853072 
EMU=2.E-1*TPI 
NT=NR*NZ 
DO  6  1=1, NT 
6  STO( I )=0. 

DO  7  J=1 ,NR 
DO  7  L=1 ,NZ 
JL=(L-1 )*NR+J 
DO  8  1=1, NCT 

IF ( XNTC ( I )  .NE.  1.)  GO  TO  14 
DZC=0. 

GO  TO  15 

14  DZC  =  (ZCR(I)-ZCL(I))  /  (XNTC(I)-l) 

15  NTCPAR  =  XNTC(I)  +  .0001 

PREVIOUS  STATEMENT  CHANGED  11-23-82  TO  FLOATING  -  D.V. 

DO  9  I K = 1 ,  NTCPAR 
R2=R(J)/RCC(I ) 

Z2=-( (ZCL (I )+( IK-1 )*DZC)-Z(L ) )/RCC ( I ) 

A=CCC(I )*EMU*R(J) 

A1=CCC( I )*EMU/RCC ( I ) 

C1=SQRT( (l.+R2)**2+Z2*Z2) 

C2=(1.-R2)**2+Z2*Z2 

AK=SQRT(4.*R2)/C1 


r# 


ft 


m 


$ 


0 


107  FORMAT  ('  MESH  POINT  COINCIDENT  WITH  COIL') 
K2=AK*AK 
CK2=1.-K2 

CALL  CELI2(E,AK,1.,CK2,IER) 

CALL  CELI2(K,AK,1.,1.,IER) 

IF (R2.NE .0. )  C3=2./(TPI*AK*SQRT (R2) ) 
IF(IFIELD)  10,11,12 

10  BRI=-K+E*(1.+R2*R2+Z2*Z2)/C2 

IF (R2.NE.0. )  BRI=A1*BRI*Z2/(TPI*R2*C1 ) 

IF (R2.EQ.0. )  BRI=0. 

STO(JL )=STO(JL )+BRI 
GO  TO  13 

11  BZI=K+E*(1.-R2*R2-Z2*Z2)/C2 
BZI=A1*BZI/(TPI*C1) 

STO( JL )=STO( JL )+BZI 
GO  TO  13 

12  IF(R2.NE.O.)  PSI=A*C3*((1.-.5*K2)*K-E) 

IF (R2.EQ.0. )  PSI=0. 

ST0( JL )=ST0(JL )+PSI 

13  CONTINUE 
9  CONTINUE 
8  CONTINUE 
7  CONTINUE 

RETURN 

END 


SUBROUTINE  CELI2(RES,AK,A,B,IER) 
I ER =0 
CK=AK*AK 
IF(CK-1.)  7,7,5 
5  IER-1 
GO  TO  100 
7  GE0=SQRT ( 1 ,-CK ) 

IF(GEO)  50,10,50 
10  I F ( B )  20,40,30 
20  RES=-1 . E 3 7 
GO  TO  100 
30  RES=1.E37 
GO  TO  100 
40  RES=A 
GO  TO  100 
50  AR 1=1. 

AA=A 

AN=A+B 

W=B 

60  W=W+AA*GEO 
W=W+W 
AA=AN 
AAR  I =AR I 


0 


IF  <  AAR  I -GEO-1 .E-4*AARI )  80,80,70 
70  GE0=SQRT(GE0*AARl) 

GE0=GE0+GE0 
GO  TO  60 

80  RES=.78539816*AN/ARI 
100  RETURN 
END 


SUBROUTINE  ECORR(STO) 

COMMON/A/  RB(50) ,ZB ( 50) ,RF(50) ,ZF ( 50) 

DIMENSION  ST0(1 ) ,CM(100) ,PN( 100) ,NRF (50) ,MZF (50) ,EH0( 100) 

COMMON/B/RCC ( 100 ) , ZCL ( 100 ) , ZCR ( 1 00 ) ,  XNTC ( 100 ) , CCC ( 100 ) 

COMMON/C/R(200),Z(200),NR,NZ,NCT 

DIMENSION  BE ( 100), EHNM( 1000), IN(100),ARR( 1000) 

DIMENSION  BZC( 100) ,BZA(100) 

REAL  K,K2 

TPI =6. 283185072 

EMU=2.E-1*TPI 

READ(2,*)  ME 

DO  40  1=1, ME 

READ(2,*)  NRF ( I ) ,MZF ( I ) 

40  CONTINUE 

DO  41  1=1, ME 
READ(2,*)  EHO(I) 

41  CONTINUE 
READ(2,*)  NE 
READ(2,*)CM0 
READ(2,*)PN0 
DO  42  1=1, NE 
READ(2,*)  RB ( I ) , ZB ( I ) 

42  CONTINUE 

DO  24  M=1,ME 
MM=MZF(M) 

NM=NRF(M) 

RF (M)=R(NM) 

24  ZF(M)=Z(MM) 

PRINT  43,CM0,PN0 

43  FORMAT (// 1  CMO= ' , 1PE13.5, '  PNO= ' , 1PE13.5) 

DO  25  M=1,ME 

25  CM(M)=CMO 

DO  26  N=1 ,NE 

26  PN(N)=PNO 
PRINT  110 

110  FORMAT ( // '  FIELD  CORRECTION  PARAMETERS’) 

PRINT  111 ,NE ,ME 

111  FORMAT(//, '  NE  = 1 , 1 12 , '  ME = ' , 1 2 , ’  EHO(M)  = ' , 1PE10.3) 
PRINT  112 

112  FORMAT (//, 1  RADIAL  MESHS  FOR  FIELD  CONTROL',/) 

PRINT  113, (RF ( I ) , 1=1 ,ME ) 


114  FORMAT (//,'  AXIAL  MESHS  FOR  FIELD  CONTROL',/) 

PRINT  113, (ZF ( I ) , 1=1, ME) 

PRINT  115 

115  F0RMAT(//10X, '  CURRENT  LOOPS ' ,/,10X, 'RADII ' ,/) 
PRINT  113, (RB( I ) , I=1,NE) 

PRINT  116 

116  F0RMAT(/,10X, 'AXIAL  POSITION',/) 

PRINT  113, ( ZB ( I ) , 1=1, NE) 

DO  27  N=1,NE 
DO  27  M=1,ME 
R2=RF (M)/RB(N) 

Z2=-(ZB(N)-ZF (M) )/RB(N) 

A1=EMU/RB(N) 

C1=SQRT( (l.+R2)**2+Z2*Z2) 

C2=(1.-R2)**2+Z2*Z2 
AK=SQRT(4.*R2)/C1 
IF(C2.EQ.O.)  PRINT  119 

119  FORMAT ( '  GRID  POINT  COINCIDENT  WITH  CURRENT  LOOP') 
IF(C2.EQ.O.)  CALL  EXIT 
K2=AK*AK 
CK2=1.-K2 

CALL  CELI2(E,AK, 1. ,CK2, IER) 

CALL  CELI2(K,AK,1. ,1. ,IER) 

NM=(M-1)*NE+N 

EHNM(NM)=K+E*(1.-R2*R2-Z2*Z2)/C2 

27  EHNM(NM) =EHNM(NM)*A1/ (TPI*C1) 

DO  28  N=1,NE 

BE(N)=0* 

DO  29  M=1 ,ME 
I 1=NRF (M) 

J1=MZF (M) 

I1J1=(J1-1)*NR+I 1 
NM=(M-1)*NE+N 

29  BE(N)=BE(N)+(EH0(M)-ST0(I1J1) )*EHNM(NM)*CM(M) 

28  CONTINUE 

DO  30  1=1, NE 
DO  30  J=1 ,NE 
IJ=(J-1)*NE+I 
ARR(IJ)=0. 

DO  31  M=1,ME 

IM=(M-1)*NE+I 

JM=(M-1)*NE+J 

31  ARR(IJ)=ARR(IJ)+EHNM{IM)*EHNM<JM)*CM(M) 

30  CONTINUE 

DO  32  1=1, NE 
1 1  =  ( I-1)*NE+I 

32  ARR( 1 1 ) =ARR( 1 1 )+PN( I )*RB( I ) 

CALL  SAXB(NE,NE,1,ARR,BE,0, IN.KER) 


IF(KER.NE.O.)  CALL  EXIT 
PRINT  121 

121  FORMAT (////, '  CORRECTION  CURRENTS',/) 

PRINT  113, (BE(N) ,N=1,NE) 

DO  33  M=1 ,ME 
I 1=NRF (M) 

J1=MZF (M) 

I1J1  =  (J1-1)*NR+1 1 

33  BZC(M)=EHO(M) -STO( I 1 J 1 ) 

DO  34  M=1,ME 
BZA(M)=0. 

DO  34  N=1,NE 
NM=(M-1)*NE+N 

34  BZA(M)=BZA(M)+EHNM(NM)*BE(N) 

PRINT  122 

122  FORMAT(//,'  CORRECTION  FIELD  REQUESTED',/) 
PRINT  113, (BZC(M),M=1,ME) 

PRINT  123 

123  FORMAT (//, '  ACTUAL  FIELD  RETURNED*,/) 

PRINT  113, (BZA(M) ,M=1,ME) 

L=NCT 

DO  35  1=1, NE 
L=L+1 

RCC(L)=RB( I ) 

ZCL(L)=ZB( I) 

ZCR(L)=ZB(I) 

CCC(L)=BE(  I ) 

3j  XNTC(L)=1 
NCT=NCT+NE 
RETURN 
END 


e 
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The  method  of  characteristics  and  multiple-scaling  perturbation  techniques  are  used  to 
study  the  space-charge  instability  of  electron  beams.  It  is  found  that  the  stable  oscillating 
state  (virtual  cathode)  created  when  the  space-charge  limit  is  exceeded  is  similar  to  a  col¬ 
lisionless  shock  wave.  The  oscillatory  solution  originates  at  the  bifurcation  point  of  two 
unstable  steady  states.  Complementary  behavior  (virtual  anode)  results  when  an  ion  beam 
exceeds  its  space-charge  limit.  The  virtual  cathode  can  also  exist  in  the  presence  of  a  neu¬ 
tralizing  heavy-ion  background.  The  Pierce  instability,  where  the  electron  and  ion  charge 
densities  are  equal,  is  a  special  case  of  this  broader  class.  Estimates  of  the  nonlinear  growth 
rate  of  the  instability  at  the  space-charge  limit  are  given. 


I.  INTRODUCTION 

Since  the  discovery  of  the  Child-Langmuir  rela¬ 
tion l,:  it  has  been  known  that  exceeding  the  limiting 
current  of  a  diode  leads  to  the  development  of  a  vir¬ 
tual  cathode.  Subsequently,  numerous  papers  were 
written  on  experiments  and  theory  relating  to 
space-charge— limited  flows.  Reference  3  provides 
an  excellent  background  and  bibliography.  More  re¬ 
cently,  the  exact  steady-state  solutions  for  electron 
beams  in  one-dimensional  relativistic  diodes4  and 
bounded  drift  spaces5  were  derived.  It  is  easily  seen 
that  for  sufficiently  large  currents  there  exist  two 
steady  states  for  an  electron  beam,  only  one  of 
which  is  stable.*1  At  the  space-charge  limit  (SCL) 
these  two  states  coalesce,  and  above  the  SCL  they 
disappear.  As  current  is  increased  past  the  SCL, 
the  beam  develops  a  jump  instability  and  relaxes 
into  an  oscillating  state. 

In  the  early  1960’s,  computer  models  were 
developed  which  quantitatively  depicted  the  non¬ 
linear  oscillatory  nature  of  the  virtual  cathode/-'1’ 
These  were  one-dimensional,  nonrelativistic,  electro¬ 
static,  multiple-sheet  models.  References  8  and  9 
qualitatively  pointed  out  many  interesting  dependen¬ 
cies  of  the  oscillation  frequency  and  potential 
minimum  position  on  injected  current,  thermal 
spread,  and  circuit  resistance.  Reference  10  presents 
computer  experiments  with  one  and  two  species. 

The  phenomenon  of  virtual-cathode  formation  in 
intense  relativistic  electron  beams  figures  prom¬ 
inently  in  a  number  of  high-interest  research  areas. 
Devices  used  to  produce  high-current  ion  beams  for 
inertial-confinement  fusion — pinch  reflex  diodes11  12 
and  reflex  triodes1’  14 — depend  on  the  virtual 
cathode  to  inhibit  electron  transport  and  use  its  po¬ 
tential  well  to  accelerate  ions.  The  recent  concept  of 
the  spherical  electron-to-ion  converter'5  requires  a 


virtual  cathode. 

The  virtual  cathode  plays  a  dominant  role  in  areas 
other  than  production  of  light-ion  beams  for  fusion. 
It  is  attributed  with  the  main  role  in  collective  ion 
acceleration  in  neutral  gas.16,17  Control  of  virtual- 
cathode  motion  is  the  mechanism  for  acceleration  in 
the  Ionization  Front  Accelerator.18'11’  It  is  also  the 
acceleration  method  in  two  concepts  for  collective- 
effect  accelerators.20  *1  A  final  application  is  the  use 
of  virtual-cathode  oscillations  to  produce  high- 
power  short-wavelength  microwaves.22-25  Experi¬ 
ments  using  reflex  triodes  have  already  produced  1.4 
GW  of  power  with  12%  beam-to-rf  energy- 
conversion  efficiency.26 

In  this  paper  we  use  multiple-scaling  perturbation 
techniques  to  study  the  time-dependent  behavior  of 
a  beam  when  the  SCL  is  exceeded.  We  derive  esti¬ 
mates  for  the  nonlinear  growth  rate  of  the  ensuing 
instability  and  show  that  even  below  the  SCL  the 
beam  is  unstable  to  sufficiently  large  perturbations. 
The  method  can  be  applied  to  a  wide  class  of  prob¬ 
lems,  but  here  we  treat  the  short-circuited  one¬ 
dimensional  electrostatic  diode  depicted  in  Fig.  1  as 
the  simplest  model  containing  the  appropriate  phys¬ 
ics.  We  show  that,  at  least  in  one  dimension,  an  ar¬ 
bitrary  heavy-ion  background  does  not  alter  the 
qualitative  behavior  of  the  beam  and  present  numer¬ 
ical  results  that  exhibit  virtual-cathode  oscillations 
for  a  neutral  beam. 

II.  PHYSICAL  DESCRIPTION 

Simulations  were  carried  out  in  conjunction  with 
the  theory  presented  in  the  next  section  using  a 
two-dimensional,  relativistic,  electrostatic,  particle- 
in-cell  code.  The  code  can  solve  self-consistently  for 
the  time-dependent  trajectories  of  tens  of  thousands 
of  plasma  particles  over  thousands  of  plasma 
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FIG.  1.  Schematic  of  the  short-circuited  one- 
dimensional  electrostatic  diode  modeled  in  this  study. 


periods.  All  variables  are  expressed  in  dimensionless 
terms.  Therefore  length  is  in  units  of  c/wp,  time  is 
measured  in  units  of  and  particle  velocity  is 
given  by 

t'i  =  Pi Y  It  =  1.2,3), 

where  cop  is  the  initial  electron  plasma  frequency. 

In  these  simulations  a  monoenergetic  51-keV  elec¬ 
tron  beam  is  injected  into  a  Cartesian  geometry. 
The  left  and  right  boundaries  are  grounded, 
representing  a  planar  short-circuited  diode.  Periodic 
boundary  conditions  in  the  transverse  direction 
make  configuration  space  effectively  one¬ 
dimensional.  In  general,  the  simulation  had  62  cells 
in  the  longitudinal  direction  modeling  a  length  of 
1.0 c/iop.  The  time  step  was  0.0125wp-1.  Twenty 
particles  were  injected  per  cell. 

A  detailed  discussion  of  the  physical  dynamics  of 
the  virtual  cathode  based  on  these  numerical  results 
is  appropriate  here.  The  usual  graph  of  potential 
minimum  6m  in  the  diode  versus  electron-beam 
current  a  is  shown  in  Fig.  2.  The  parameter  a  will 
be  discussed  later.  When  a  is  increased  above  the 
space-charge  limit,  6„  jumps  from  the  stable 
normal-C  branch  to  the  oscillatory  stable  branch. 
The  amplitude  and  position  of  d>„  while  on  the  os¬ 
cillatory  branch  describe  a  limit  cycle,  as  expected 
for  a  relaxation  oscillation  which  this  represents. 
Typical  limit  cycles  are  depicted  in  Fig.  3.  As  a  is 
increased  further.  6m,  the  oscillation  frequency,  and 
virtual-cathode  position  within  the  diode  asymptoti¬ 
cally  approach  limiting  values.  If  a  is  decreased,  the 
oscillation  amplitude  decreases,  and  the  posi- 
,ion  of  <t>m  moves  toward  the  diode  center.  The  elec¬ 
tron  flow  reverts  to  the  equilibrium  steady  state 
when  the  perturbation  due  to  the  rate  of  change  of 
diode  current  below  the  space-charge  limit  is  suffi¬ 
ciently  large.  This  normally  occurs  before  the  bifur¬ 
cation  point  is  reached.  The  entire  process  forms  a 
hysteresis  loop,  which  is  depicted  in  Fig.  2. 

The  virtual  cathode  originates  at  the  bifurcation 
point.  This  is  the  intersection  of  the  oscillatory  state 
with  the  C-overlap  and  partially  reflected  solution 
branches.'  The  bifurcation  point  cannot  be  reached 


BIFURCATION 

POINT 


C  UVKRLAPyj 

■  ,  *: 


HYSTERESIS  LOOR 


FIG.  2.  Electrostatic  potential  minimum  as  a  function 
of  current  a.  Plot  depicts  the  various  possible  solutions 
such  as  normal-C  flow  (stable).  C-overlap  (unstable),  the 
partially  reflected  solution  (unstable),  and  the  oscillatory 
virtual  cathode  (stable).  Motion  around  the  hysteresis 
loop  is  denoted  by  arrows. 


in  the  short-circuited  diode.  Of  the  three  branches 
emanating  from  it,  two  (the  steady  ones)  are  physi¬ 
cally  unstable  while  the  oscillatory  branch  is  numer¬ 
ically  unstable  at  this  point.  This  results  because  the 
limit  cycle  at  the  bifurcation  point  is  infinitesimally 
small,  so  that  simulation  codes  lose  resolution  before 
it  can  be  reached.  Loss  of  resolution  creates  a 
small-amplitude  high-frequency  oscillation  observed 
in  this  study  and  previously.1*  This  result  is  numeri¬ 
cal,  not  physical. 

This  problem  can  be  overcome  if  we  eliminate  the 
hysteresis  loop.  Then  the  C-overlap  branch  disap¬ 
pears,  and  we  can  get  to  the  bifurcation  point  along 
the  stable  normal-C  branch.  This  can  be  accom¬ 
plished  in  several  ways.  The  most  appropriate  in 
this  study  is  to  have  a  retarding  potential  difference 
across  the  diode  equal  to  the  injected-electron  kinet¬ 
ic  energy.  Then  the  C-overlap  solution  vanishes, 
and  the  bifurcation  takes  place  at  the  space-charge 
limit  astL-  The  oscillation  can  be  described  as  a 
small  perturbation  on  the  beam  rather  than  the  radi¬ 
cal  change  that  results  in  the  short-circuited  diode 
when  aSci  *s  exceeded.  This  analysis  indicates  that 
the  onset  of  virtual -cat hode  formation  occurs  when 
the  electron  velocity  in  the  steady  state  vanishes  at 
some  point  inside  the  diode.  For  the  short-circuited 
diode,  this  occurs  at  the  diode  center;  for  the  biased 
diode,  it  occurs  at  x  </.  Because  the  oscillation  is  a 
small  perturbation  on  the  steady-state  fields  in  the 
biased  diode,  it  is  readily  seen  that  the  virtual- 
cathode  oscillation  period  at  onset  is  the  electron 
transit  time  from  the  injection  plane  to  the  position 
where  the  velocity  vanishes. 

Finally,  consider  the  particle  dynamics  during  the 
oscillation  period  for  counterclockwise  limit  cycles 
|(7<«S(  |  )  as  in  Fig.  3(c).  At  the  point  where  the 
virtual-cathode  position  is  a  minimum  and  the  po¬ 
tential  well  is  starting  to  move  to  the  right,  its  am¬ 
plitude  is  too  small  to  stop  the  electrons.  When  its 
motion  is  to  'he  left  it  opposes  the  electron  beam 
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x  (c  /  wp ) 


X  (c/Wp) 


being  responsible  for  the  onset  of  oscillatory 
behavior  in  beam  characteristics.  As  the  potential 
minimum  reaches  the  end  of  its  left  motion  the  two 
"lips”  of  the  back  reflected  stream  close.  At  this 
point  no  more  electrons  are  reflected,  and  the  well 
moves  to  the  right,  repeating  the  cycle. 

For  larger  values  of  a  the  limit  cycle  is  distorted 
into  a  figure  “8"  with  one  lobe  having  a  clockwise 
motion  and  the  other  a  counterclockwise  one  [Fig. 
3(b)].  This  transition  continues  until  the  motion  is 
completely  clockwise  [Fig.  3(a)].  It  indicates  a 
change  in  the  particle  bunching  process  and  is  relat¬ 
ed  to  the  fraction  of  current  which  is  reflected 
versus  transmitted  from  the  injected-electron  beam. 
In  Fig.  3(a)  most  of  the  beam  is  reflected,  whereas  in 
Fig-  3(c)  most  electrons  are  transmitted. 

III.  THEORY 

The  one-dimensional  motion  of  electrons  in  the 
diode  is  given  by  the  equations  of  continuity  and 
momentum  conservation  for  the  electrons  plus 
Poisson's  equation.  They  are  expressed  here  as 

p-+(pv)~  —  0,  (la) 


FIG  3.  Typical  virtual-cathode  limit  cycles  in  the  clas¬ 
sical  short  circuited  one-dimensional  diode  with  an  injec¬ 
tion  energy  of  51  keV.  *a>  «  2.5,  'b'  «  =  2.0.  and  icl 

a  1.4.  Motion  in  •  a >  is  clockwise  and  in  ic)  is  counter¬ 
clockwise  / 

and  causes  particle  bunching,  Since  the  well  is 
deeper,  the  stream  velocity  will  vanish  at  some  loca¬ 
tion  and  then  become  negative.  Here,  the  second 
derivative  of  the  velocity  td:u  /dx:)  is  also  negative. 
In  this  process  the  stream  is  continuously  deformed 
to  create  a  double-valued  negative  velocity  pro¬ 
trusion.  The  entire  system  is  three-valued  'Fig.  4), 
as  in  a  eollisionless  shock  wave.-’  Here,  the  region 
of  triple  flow  is  not  limited,  as  in  usual  collisionless 
shocks,  by  the  presence  of  a  transverse  magnetic 
field, but  rather  by  the  presence  of  the  walls. 
Indeed,  the  reflected  part  detaches  from  the  main 
beam  and  exits  through  the  anode  periodically,  thus 


v7+vvx  =  ~  ie/m)£  ,  (lb) 

e(£x=P+Ph  .  dc) 

where  0<x<I.  The  subscripts  denote  differentia¬ 
tion  with  respect  to  that  variable.  Electron  and 
heavy-ion  charge  densities,  and  electron  velocity, 
electric  field,  time,  and  position,  are  indicated  by  p, 
ph,  v,  c,  /.  and  Jc,  respectively.  The  appropriate 
boundary  conditions  are  i>(0,f)  =  r0,  p(0,t)—p0,  and 
J  £dx=  0.  An  ion  component  is  placed  in 
Poisson's  equation  in  order  to  discuss  two-species 
space-charge  flow.  Conservation  equations  of  mass 
and  momentum  for  the  heavy  ions  are  not  included, 
because  it  is  assumed  that  their  velocity  does  not 
change  appreciably  during  their  passage  through  the 
diode. 

In  order  to  simplify  the  mathematics,  we  intro¬ 
duce  dimensionless  equations  for  our  model.  They 
are 

n,  < nu  I,  -0  ,  (2a) 

u,  mii,  —  -  ctE  ,  (2b) 

Ex  =-«(  n  ■*  /i  ,  (2c  i 

w  here  n  -  p/  pn  .  u  r  /r„, 

E  -  ice, ,/m  f>„  )'  /r,,  , 

t  r,,l  //  ,  ,v  -  .v  //  , 

It  -  If  />, ,  /€„/>!  I1  ’■//!•„ 
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FIG  4  Successive  snapshots  of  electron-beam  momentum  space  and  corresponding  potential  shape  in  the  diode  for 
at  --  2.0. 1  -  1.0c  /i Time  between  frames  is  0.5 to.  Initial  beam  kinetic  energy  is  51  kev. 


fa  dimensionless  parameter  related  to  current),  and 
the  ratio  of  heavy-ion  to  electron  charge  densities, 
l-ph/  pt)  .  Alternatively,  a  may  be  written  as 
liop/vu,  where  iop  is  the  beam  plasma  frequency. 
The  boundary  conditions  for  electrons  become 
u  tO,/ i=l.  « tO,/)  =  —  1,  and  J()£cO c  =0.  For  an 
ion  beam,  /t<0,/l  =  l;  otherwise,  the  following 
derivation  is  the  same. 

A.  Unneutralized  beams 

For  an  unneutralized  beam,  setting  /  =0  and  solv¬ 
ing  by  the  method  of  characteristics2''  we  find 

n  1  =  -  t«:i /  —  v i: -t-a£()is )( r  —  vt—  1  ,  (3) 

where  v  is  the  entry  time  for  the  particle  occupying 


position  x  at  time  /,  and  £0(s)  is  the  electric  field  at 
x  =0.  The  particle  trajectories  are  found  by  utiliz¬ 
ing  Eq.  (2a),  from  which  it  follows  that 


This  yields 

x  =  ja:(/  -s)}  +  a  f  £0(s)(/  -s)ds  +  (/  —s) 


Integration  of  the  trajectory  equation  is  hard  for 
general  time-dependent  situations,  because  imposing 
the  proper  boundary  conditions  leads  to  a  nonlinear 
integral  equation  for  £,)(/)■  However,  several  special 
cases  can  be  solved  exactly.  The  problem  of  injec- 
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tion  into  an  empty  diode  can  be  integrated  until  the 
formation  of  a  singularity  in  n,  indicating  the  cross¬ 
ing  of  trajectories.'"  In  this  case,  the  stream  velocity 
becomes  three-valued  and  one  must  use  a  Vlasov- 
equation  description,'1  rather  than  system  12)  that  is 
derived  assuming  a  single  stream  of  monoenergetic 
particles.  As  described  in  Sec.  II,  this  multistream¬ 
ing  is  characteristic  of  the  oscillatory  state  created 
when  a  exceeds  its  SCL  value. 

Using  these  equations  we  can  derive  a  similar  rep¬ 
resentation  for  Fig.  2  in  terms  of  £0  and  a.  For 
steady  states,  E0it)  =  E0,  a  constant,  we  Find 

u  —  —  n  _l  =  yarU  —  s)2—aE0(t  —  s)+  1  ,  (6a) 

x  =  ja:(t  —  s)'—  ~aE0U  —s)2  +  (t  — s )  .  (6b) 

Imposing  the  conditions  x  =  1,  u  =  1  at  I  —  s  =J0, 
the  particle  transit  time,  we  note  that  /0  must  satisfy 

yror't,)  —  4-  1  =  0  .  (7) 

This  equation  has  two  positive  solutions  for 
0 <a<  coalescing  at  a  =  The  largest  one,  for 
0<cr<2v  2/3  does  not  correspond  to  a  real  flow. 
In  Fig.  5  we  show  £0  ( =  ar0/2)  vs  a.  This  represen¬ 
tation  will  be  used  in  the  discussion  of  nonlinear  sta¬ 
bility. 

B.  Neutralized  beams 

The  Pierce  instability  occurs  when  electron  and 
ion  space-charge  flow  is  considered  in  finite 
geometries  where  there  is  no  potential  difference 
across  the  boundaries.’2  The  ions  can  be  stationary 
or  moving  with  respect  to  the  electrons.  Charge 
neutrality  is  maintained  at  the  injection  plane.  This 
instability  may  have  ramifications  for  charged- 
particle  inertial  confinement  fusion  because  of  its  ef- 


OSCILLATORY 
VIRTUAL  CATHODE 


FIG.  5  Electric  field  at  the  injection  plane  vs  current  a 
for  I  0.  Plot  depicts  the  normal  C'-flovv  stable1  and  C- 
overlap  unstable'  solutions.  Oscillating  virtual-cathode 
stable1  solution  is  also  shown  Regions  I  and  II  define 
tiie  domains  of  attraction  of  (lie  norntal-C  and  virtual- 
cathode  solutions  near  the  SCL.  << 


feet  on  neuirah/ed-beam  propagation  in  the  reac¬ 
tor.1'  In  this  section,  we  show  that  the  Pierce  insta¬ 
bility  is  a  special  case  of  electron  and  ton  space- 
charge  flow.  In  general,  two-species  flow  has  steady 
and  oscillatory  states  analogous  to  one-species 
space-charge  flow  . 

The  steady-state  behavior  for  the  case  of  arbitrary 
7  can  be  found  in  a  manner  similar  to  /  =0.  Rewrit¬ 
ing  system  (2)  in  characteristic  coordinates,  we  ar¬ 
rive  at 

1-f-  +a2l  1  +  -  =0.  (8) 

n  n 

For  positive  ions  (/  >  0)  the  solution  of  (8)  after 
satisfying  the  boundary  conditions  is 

1  +1  /n  =  ( 1  +I)cos[a^/ 1  (/  —  s)] 

+  E0\/I  sinfav'/  (t  —  s)]  .  (9) 

Imposing  conditions  x  =  1,  u  =  1  at  t  —  s  =t0,  we 
find  the  system 

1=7*0-  “ITT  sin(av7/„> 

I  al  •'* 

E() 

+  — r[cos(orv//f0)— 1]  ,  (10a) 

al 

1  =  y  —  — y—  cosfaC/fo ) 


£0 

—  —  sin(av7/0) 


(10b) 


For  7  =  1,  which  implies  charge  neutralization, 
these  equations  reduce  to  the  relations  given  in  Ref. 
34  for  the  Pierce  instability.  However,  by  varying  7 
the  curves  shown  in  Fig.  6  are  obtained.  These  are 
cuts  at  constant  7  through  a  three-dimensional  con¬ 
tiguous  surface.  The  space  is  defined  by  the  axes 
£  =  £(|7’/2,  .4  =o71/:,  and  7.  The  surface  is  2 ir 
periodic  in  A  with  the  vertical  plane  at  A  =2n  being 
common  for  all  values  of  7.  For  given  A,  a  linear¬ 
ized  analysis  establishes  that  the  equilibria  denoted 
by  the  curves  are  stable  (unstable)  for  the  lowest 
i highest)  value  of  £.  At  7  =  I,  exchange  of  stability 
takes  place  at  odd  multiples  of  For  7<1,  ex¬ 
change  of  stability  occurs  at  the  points  where 
</£„/</  a  +  s. . 

It  is  evident  from  Fig.  6  that,  for  7  <  I,  there  are 
no  stable  equilibrium  solutions  in  the  neighborhood 
of  A  =  -.  Therefore  one  expects  a  virtual  cathode  to 
form  when  7  <  1  and  A  adiabat ical ly  increases  to  rr. 
We  have  found,  by  using  numerical  simulation,  that 
in  this  case  the  beam  settles  to  an  oscillatory  state, 
similar  to  the  virtual  cathode  for  unneutralized 
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FIG.  6.  Curves  of  sealed  electric  field  at  the  injection 
plane  E  vs  scaled  current  .4  for  various  values  of  charge 
neutralization  /.  Curves  represent  slices  through  a  three- 
dimensional  surface.  S  and  U  indicate  stable  and  unstable 
branches,  respectively,  for  the  /  =  1  slice. 

beams."  By  slowly  increasing  /  past  the  neutral- 
beam  value  of  1  in  our  simulation,  we  have  estab¬ 
lished  that  this  oscillation  persists.  Indeed,  finding 
this  oscillatory  state  for  /  >  1  by  other  means  would 
have  been  difficult,  because  the  simulation  would 
tend  to  follow  the  stable  steady-state  branch  that  is 
present  for  all  values  of  current. 


stable,  while  at  the  SCL  (e  =  0)  we  have  neutral  sta¬ 
bility. 

Above  the  value  a  =-,  linearized  theory  is  not 
applicable.  Utilizing  multiple-scaling  perturbation 
theory,'5  we  can  carry  out  a  nonlinear  stability 
analysis  neara=y.  In  system  (2)  we  set  n=4+f2. 
In  this  neighborhood,  perturbations  evolve  on  a 
"slow"  time  scale,  depicted  by  x=et. 

Eliminating  the  electric  field  by  combining  (2b) 
and  (2c)  and  utilizing  r,  system  (2)  becomes 


enT  +  (nu)x  =0  , 

(15a) 

(euT  +  uux  )x  =  —  ( 4  r  e1)2/! 

(15b) 

with  conditions  u(0,x)=l,  «(0,x)  = 

—  1,  and 

f  E  dx  =0  rewritten  as 

J  0 

e  J*  uTdx  +  t[u2(  1 , x )  —  tv 2( 0, x)]  =0  . 

(16) 

Substituting  the  asymptotic  expansions 

u  -  2  eiui(x,r)  +  0{ej  +  ' )  , 

(=0 

n  -  £  e'tt,(x,x)  +  0(£''  +  l)  , 

i=0 

(17) 

C.  Nonlinear  stability  analysis 

For  /<  1,  it  is  of  interest  to  establish  the  proper¬ 
ties  of  the  beam  instability  at  the  SCL — generalized 
for  I  ^  0  to  mean  the  point  where  da/dEt)  =  0.  We 
shall  carry  out  the  analysis  for  7=0,  but  our 
method  can  be  applied  to  any  similar  jump 
phenomenon. 

A  linear  stability  analysis6  about  the  steady  state 
described  by  (7 )  results  in  the  dispersion  relation 


for  u  and  n  into  (15)  and  equating  coefficients  of 
various  powers  of  e,  there  results  a  hierachy  of  equa¬ 
tions  for  the  Uj  and  nt. 

Solving  the  0(1)  system  gives 

(«0-y  )(u„+l)2  =  2(2* -l)2  ,  (18a) 

nu=  —  1  /«(>  .  (18b) 

To  solve  the  Ole)  system  we  introduce  a  new  vari¬ 
able  q  by 


i2  +  /3)e-0=2-/?-+-/?7(a2/o>  .  (11) 

where  /3=iojt().  We  have  written  the  expression  de¬ 
rived  in  Ref.  6  in  terms  of  our  dimensionless  vari¬ 
ables.  For  a  near  the  SCL  value  we  let 

a  =  4  —  ,  e  «  1  (12) 


so  that 


£1 

6 


+  1  . 


(19) 


1 20) 


and  find  from  (7)  that  near  this  value,  t0  is  approxi¬ 
mately 

fn-y-e — -  +Oirl  ,  (13) 

2  2v  2 

where  the  —  (-*-)  sign  corresponds  to  the  lower 
(upper)  branch  in  Fig.  5. 

By  substituting  in  ill),  and  assuming  (1  small,  we 
find  that 

icj  =0"  :  2v  2f  -*-Oif:)  .  1 14) 

Since  the  linearized  analysis  led  to  time  factors  of 
the  form  e"JI  in  the  perturbations,  it  follows  that  the 
lower  branch  in  Fig.  5  is  stable  and  the  upper  un- 


We  then  find 


«f> 

<7l<7  -4  ) 

U)  =C- - -  (21b) 

«(. 

with  C  a  constant  of  integration  which  is,  in  general, 
a  function  of  the  slow  time  x.  To  find  C  =C(x) 
which  determines  the  slow  evolution  of  the  pertur¬ 
bation  «|(.t,xi  we  need  to  go  to  the  next  order. 
Off2).  By  substituting  in  the  expressions  for 
un,u  |  i  and  eliminating  //  :,  we  find  that 
«:l.x, x)  satisfies 


where  r,,  is  a  constant  of  integration.  In  general, 
small  initial  perturbations  will  lead  to  the  solution 

u  Lx.t  i  -  -t-f  ^  C,(  r)ea>''ulix)  +  Ofe:)  ,  (28) 

i  -- 1 

where  <■>,  are  the  various  distinct  solutions  of  the 
dispersion  relation  U  1>  at  «=  It  is  straight¬ 

forward  to  show  that  all  modes  are  such  that 
Rev ■>,  0  except  one  for  which  w  =  0.  Thus  all  other 

modes  will  decay  in  the  fast  time  scale  and  only  the 
neutral  mode  iw=0)  will  persist.  Our  solution  after 
a  short  time  will  look  like 

q{q  -  t  ) 

u  —  Mi,  +-  € C i  r )  -  - +0(e")  .  (29) 

from  the  given  initial  conditions  it  is  easy  to 
determine  the  initial  condition  for  the  neutral  mode. 
Below  the  SCI.,  if  the  initial  conditions  are  such 
that  Ci()>  .  ic/hi,/:,  the  solution  will  evolve  to 


the  stable  lower  branch  in  Fig.  5  (region  I),  while  if 
C(0)<  —  (c/b)W2,  C— * x  in  Finite  time  (Fig.  5,  re¬ 
gion  II).  Blowup  in  finite  time  also  occurs  above  the 
SCL  for  any  C(0).  This  does  not  mean  that  the  ac¬ 
tual  solution  blows  up,  just  that  it  evolves  to  a  final 
state  far  away  from  the  two  steady-state  branches 
shown  in  Fig.  5  and  thus  is  not  accessible  by  pertur¬ 
bation  theory. 

As  can  be  seen  in  (25),  the  blowup  above  the  SCL 
is  described  by  a  tangent  function,  therefore  the 
growth  rate  we  find  for  this  case  must  be  appropri¬ 
ately  interpreted.  Note  that  the  linearized  dispersion 
relation  seems  to  suggest  an  imaginary  exponential 
growth  rate  above  the  space-charge  limit.’  6  In  view 
of  our  results,  we  see  that  this  is  actually  misleading. 
Moreover,  we  find  that  even  below  the  SCL  the 
stable  steady-state  branch  can  be  destabilized  by  suf¬ 
ficiently  large  perturbations. 

Our  results  agree  with  the  linear  theory,  provided 
we  consider  the  limit  where  the  latter  becomes  ap- 
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plicable.  Thus  we  must  compare  the  linear  theory 
with  (26)  as  r—+  x  (near  the  stable  branch)  and 
(27)  as  r— *—  x  (near  the  unstable  branch).  To 
demonstrate  this  we  set  r  =  et  in  (26)  and  consider 
the  limit  r— »  +  x  .  Then 

Cil)=~~ tanh  V  C^  (ft  -f  r0) 
o  a 


1  —  exp 

2v  cb  ,  , 

a 

2  v  cb  . 

1  +  exp 

-  i€t  +  r0) 

a 

Y  1—2  exp  — -^-^-(et+r0) 
b  a 


+  ■  ■  ■ 


(30) 


We  see  that  the  decay  rate  of  the  perturbations  as 
/  — ►  +  x  is  equal  to 


which  is  the  same  as  that  found  by  the  linearized 
analysis.6  This  quantity  is  important  as  it  also 
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determines  an  initial  “growth  rale”  for  the  jump  in¬ 
stability  above  the  SCL  described  by  (25).  If  we  sub¬ 
stitute  the  original  dimensional  time  variable  into 
our  expressions  and  write  the  deviation  of  a  from  its 
value  at  the  SCL  as 

e  =  (a-ascl  )l/’  =  (a— 4  )l/2  .  (32) 

we  find  that  the  "growth  rate"  is  given  by 


a  / 


i  1/2 

t'n  3  / 


(33) 


(34) 


Of  course,  for  the  expansions  in  (28)  to  be  valid,  we 
must  have  C(r) «  1/e.  However,  while  C  is  not  too 
large,  (25)  gives  a  reliable  estimate  for  the  growth 
rate  of  the  instability. 
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Abstract 


TABLE  1 


High  power  microwaves,  up  to  Gigawatt  levels  in 
the  centimeter  regime,  have  been  observed  in  reflex 
triode,  foil  and  foil  less  diode  systems.  Generation 
efficiencies  range  from  1%  to  121.  The  source  of  the 
microwaves  is  an  oscillating  virtual  cathode  -  the 
nonlinear  state  which  develops  when  the  electron  beam 
injection  current  exceeds  the  space-charge  limiting 
cj^rent  defined  by  the  beam  energy  and  wave  guide  geo¬ 
metry.  This  stable  oscillation  results  in  severe 
longitudinal  charge  bunching  giving  rise  to  large  time 
dependent  current  variations.  The  experimental  fre¬ 
quency  dependence  and  broadband  characteristics  are 
explained  by  the  scaling  of  the  oscillator  frequency 
w'tn  vnc  'y  ,  wnere  nb  is  the  beam  density  and  y  its 
relativistic  factor,  in  conjunction  with  the  Oiild- 
^angmu'r  relation.  The  optimal  design  for  a  narrow- 
band  millimeter  wave  vircator  is  based  on  a  foilless 
d'ode  with  a  strong  axial  magnetic  field.  It  will  be 
tunable  over  an  order  of  magnitude  in  frequency  by 
varying  the  magnetic  field  strength. 

Introduction 


HIGH  POWER  MICEOWAVE  OENEEATION  HAJ  HEN  WTTNIMED  IN  VAEIOUI 
ELECTRON  BEAM  CONFIGURATIONS  WHEN  VIRTUAL  CATHODE!  ARE  FORMED 
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asymptotical ly  approach  a  limiting  value  for  current 
above  the  space-charge  limit.  Second,  the  fundamenta' 
oscillation  frequ-  approximately  the  relativ¬ 
istic  beam  plasr  __  given  by 


0*  the  several  millimeter  sources  that  are  in 
various  stages  of  development,  the  virtual  cathode 
oscil’ator  (Vircator)  nas  a  combination  of  character¬ 
istics  wmch  recommend  it  for  high  frequency  use. 
First,  tne  frequency  of  the  vircator  is  tunable  by 
cnanging  the  magnitude  of  an  imposed  axial  magnetic 
field,  eliminating  any  requirement  to  change  the  phy¬ 
sical  structure  of  tne  device.  A  single  vircator  will 
be  tunable  over  an  order  of  magnitude  in  frequency 
(e.g.,  10  Ghz  -  100  Ghz).  Second,  tne  bandwidth  of 
tne  generator  can  be  narrow  or  broad  based  on  magnetic 
field  shaping  and  the  use  of  beam  limiters  described 
below.  Third,  because  the  vircator  functions  above 
the  space-charge  limiting  current  for  the  electron 
bean,  give"  efficient  operation,  it  should  be  capable 
of  much  higher  power  tnan  other  microwave  sources. 
Finally,  tne  lack  of  passive  resonating  structures  to 
produce  the  transmitted  wave  reduces  the  problem  of 
*ieiq  e-Tssion.  This  also  increases  the  maxi-inn  pos¬ 
sible  generator  power. 

Experimentally,  the  virtual  cathode  has  already 
p-Oven  itself  to  be  a  copious  microwave  source.1* 
Tab’e  1  lists  experiments  which  have  been  carried  out 
to  late.  With  the  exception  of  the  Didenko  experiment 
at  *omsx,"  the  frequency  spectra  have  all  had  a  broad 
bandwidth  and  relatively  low  efficiency.  Neverthe- 
'ess,  even  at  low  efficiency  the  experiment  at  Harry 
Diamond  Laborat jr-ies5  using  a  foilless  diode  (1  MV,  30 
•A;  produced  Gigawatts  of  power  in  the  Ku  band.  It  is 
one  of  the  most  powerful  centimeter  wavelength  micro- 
wave  sources  available.  It  will  be  shown  later  that 
tne  foil  less  diode  in  a  shaped  axial  magnetic  field 
represents  the  best  configuration  for  a  high  frequency 
dev i ce. 


Microwave  Generation 


Although  only  scaling  relations  are  presently 
avsi'aole,  qualitative  dependencies  of  virtual  cathode 
parameters  on  bean  kinetic  energy  and  injected  current 
are  known.9.9  First,  potential  amplitude,  position, 
and  oscillation  frequency  all  have  the  same  functlona’ 
dependence  on  injected  beam  current.  These  parameters 
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m.  e 
b 
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,1/2 


(1) 


where  n°  is  the  electron  beam  number  density  at  injec¬ 
tion  yq  is  the  beam  relativistic  factor,  e  is  the 
electron  charge  and  m  is  its  mass.  In  particular, 
the  oscillation  frequency  from  one-dimensional  elec¬ 
trostatic  and  two-dimensional  electromagnetic  nme-i- 
cal  simulations  varies  such  that 
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The  value  of  /ZiT  is  an  empirical  result  which  has  not 
yet  been  derived  theoretically.  The  value  of  uos; 
increases  with  current  monotoni cal ly .  If  the  injec¬ 
tion  current  exceeds  the  space-charge  limiting  current 
by  a  factor  of  three  or  greater,  u)0SC  is  close  tu 
the  maximum  value.  Equation  11  in  conjunction  with 
the  Child-Langmuir  law  describing  space-charge  limited 
diode  emission  explains  the  experimental  linear  depen¬ 
dence  of  frequency  on  the  square  root  of  diode  voltage 
in  foil  diodes  and  reflex  triodes. 


The  virtual  cathode  oscillates  stably  at  a  set 
frequency  in  both  time  and  space.  This  fluctuating 
potential  barrier  acts  as  a  gate  to  reflect  some  elec¬ 
trons  and  transmit  others.  The  motion  of  the  gate 
bunches  charge.  In  two  dimensions  the  charge  bunch 
and  virtual  cathode  are  separated  spatially. 

By  analogy  this  configuration  represents  an  LC 
oscillator.  The  virtual  cathode  acts  as  a  capacito- 
to  store  the  beam  kinetic  energy.  During  that  portior 
of  the  limit  cycle  in  which  the  potential  is  greater 
tnan  the  injected  bean  energy,  charge  is  constrained 
to  remain  near  the  anode.  This  starves  the  virtual 
cathode  so  that  its  amplitude  decreases  below  (y0 
-1)  mc2/e.  Once  this  occurs  the  charge  bunch  is 
transmitted.  The  electron  motion  represents  a  large 
time  varying  current  through  an  inductor.  The  pres¬ 
ence  of  charge  away  from  the  anode  reestablishes  the 
virtual  cathode,  and  the  cycle  repeats.  The  effect  on 
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beam  current  can  be  examined  using  simulations.  A  net 
current  diagnostic  is  given  In  Fig.  1  where  the  probe 
is  positioned  between  the  anode  and  virtual  cathode. 
The  injected  current  in  units  of  mc3/e,  v0,  Is  3.4 
times  the  limiting  current,  vt,  and  Yo  15  3.5. 
Note  that  the  virtual  cathode  can  actually  reverse  the 
direction  of  current.  The  average  current  value  is 

Vf 

The  oscillating  current  generates  micro- 
waves.10’11  The  wave  frequency  Is  the  oscillation 
frequency  of  the  virtual  cathode.  The  wave  propagates 
down  the  drift  tube  in  i  W  waveguide  mode,  which 
determines  the  wavelength  and  phase  velocity  of  the 
wave  in  the  guide.  The  field  configuration  Is  evi¬ 
dent  in  simulations  where  there  is  no  axial  magnetic 
field.  If  a  cold  beam  is  injected  and  azimuthal  sym¬ 
metry  is  assumed  by  tne  code,  the  only  nonzero  fields 
are  E,,  Er,  and  Be  in  cylindrical  geometry. 
These  three  fields  define  a  TM  wave  traveling  in  the  z 
direction. 


Figure  1.  he:  currert  amp’-tude  ard  sDectrur  fen  a 
probe  placed  between  the  anode  and  v-rtua1  cathode. 
•0=3.5,  .c*3.4.,  • 


The  largest  impediment  to  constructing  an  effi¬ 
cient  vircator  is  the  effect  which  heating  has  on  the 
microwave  generation  efficiency.  As  noted  earlier12, 
beam  temperature  significantly  damps  out  the  amplitude 
of  the  potential  oscillation.  This  can  be  understood 
in  the  following  way.  For  a  nonoenergetic  beam  all  of 
the  charged  particles  bunch  at  the  sane  location. 
Mathematically  this  represents  a  singularity  where  the 
charge  density  goes  to  infinity.  In  reality  the 
charge  bunch  is  not  infinitely  dense,  but  it  does 
become  several  times  greater  than  the  initial  bean 
injection  density.  The  severity  of  the  charge  bunch¬ 
ing  leads  to  efficient  microwave  generation.  If,  on 
the  other  hand,  the  beam  has  a  spread  in  axial  momen¬ 
tum,  the  electrons  will  stop  at  different  locations  in 
the  potential  well.  This  tends  to  limit  the  charge 
bunching  and  the  amplitude  of  the  osc i 1 1  at i ng  electric 
and  magnetic  fields.  The  effect  of  beam  temperature 
in  reducing  the  RF  efficiency  of  the  vircator  has  been 
witnessed  in  one-dimensional  electromagnetic  simula¬ 
tions.13  A  beam  spread  of  less  than  3i  in  energy 
reduces  the  microwave  generation  efficiency  from  2W 
to  approximately  21.  Under  these  conditions  the 
vircator  is  nothing  more  than  a  Barkhausen  oscil¬ 
lator.1*  Indeed,  the  low  efficiency  and  broad  band¬ 
width  observed  In  most  of  the  experiments  to  date  can 
probably  be  attributed  to  the  effects  of  electron  re- 
flexing  in  the  diode  region  resulting  in  beam  heating. 


Vircator  Design 


Based  on  ou'  current  theoretical  knowledge  o‘  tne 
virtual  cathode,  numerical  simulations,  and  experimen¬ 
tal  results,  the  following  points  must  be  considered 
in  the  design  of  a  coherent,  high  frequency  vircator. 
First,  generation  of  high  microwave  'reguencies  re¬ 
quires  1 arge  beam  densities.  Extremely  nigh  bea- 
densities  ( nP  >  101*  cm-3)  have  been  obtained  from  a 
foil  less  diode.  The  beam  plasma  frequency  scales  lin¬ 
early  with  the  electron  cyclotron  frequency  due  to  tne 
magnetic  field.15  This  Is  of  significance,  because  it 
Implies  that  a  single  vircator  can  be  tuned  over  an 
order  of  magnitude  in  frequency  (eg.  10-100  GHz)  sim¬ 
ply  by  changing  the  axial  magnetic  field  strength 
without  changing  the  physical  structure  of  the  device. 

Second,  both  the  oscillation  frequency  and  net 
current  asymptotically  approach  a  value  as  injected 
current  is  increased  above  the  space-charge  limit. 
Thus,  using  a  very  large  value  of  v0/vt  aoe5  not 
substantially  increase  frequency  or  RF  efficiency 
where  efficiency  is  given  by 

n  <  — %*— .  (3 

(v0-l) 

A  foilless  diode  in  a  strong  axial  magnetic  f’e’d  pro¬ 
duces  a  very  thin  annular  beam.  Since  fpr  ar  an¬ 
nular  bean  is  larger  than  for  a  solid  beam  o*  tne  sa-e 
area,  the  value  of  v0/»£  will  be  sna'ier  for  tne 
same  beam  current. 


In  order  to  assjre  narrow  bauw'dC',  nig- 
efficiency  microwave  generation  at  higr  or  low  fre¬ 
quency  from  the  vircator,  the  following  cnara::er- 
istics  must  be  met.  First,  no  reflexing  of  ele::rp's 
in  the  region  between  the  real  and  virtual  cat-pdes 
must  occur.  Any  axial  magnetic  field  must  be  shape: 
to  divert  the  electrons,  or  flux  excluders  ~x-st  be 
employed  to  confine  the  magnetic  field  to  the  d-ode 
region.  The  latter  arrangement  will  allow  the  vya’ 
space-charge  electric  field  to  perfo-n  the  ro’e  :• 
expelling  electrons  to  the  waveguide  wa’l.  In  add-- 
tion,  if  the  beam  is  annular,  a  collimator  may  be  use: 
to  help  prevent  reflexing  of  electrons  ha:>  to  tv 
cathode.  Second,  the  electron  beam  must  be 
Experimental15  and  theoretical15  results  indicate  that 
foilless  diodes  create  low  enittance  beams.  ta"r3' 
flow,  where  the  electron  larmor  orbit  is  sna’ier  t',a' 
the  bean  thickness,  is  obtained  when16 
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where  a  is  the  orbit  radius,  6  is  the  radial  spacing 
between  the  cathode  and  drift  tube  wall  (which  acts  as 
the  anode)  and  wc  is  the  electron  cyclotron  fre¬ 
quency  given  by  eBz/mc.  Low  beam  scatter  is  also 

assured,  because  of  the  lack  of  a  foil.  Fina’ly,  the 
diode  voltage  and  injected  current  must  be  constant. 
More  appropriately  stated,  the  impedance  must  be  con¬ 
stant.  Flat-top  voltage  pulses  can  be  attained  m  a 
variety  of  ways  in  several  diode  configurations.  How¬ 
ever,  at  high  voltages  the  foilless  diode  operates  as 
a  purely  resistive  load,  therefore  «osc“  'V  =  1  >"? 
is  constant.  Also,  absence  of  diode  closure  m  some 
f oi  1  less  d’ode  experiments  makes  a  long  pu'se  de<-a 
possible. 

It  is  evident  from  this  discussion  tnat  tne 
foilless  diode  in  a  strong  axial  magnetic  field  re.- 
resents  the  optimal  configuration  for  a  higr  fre¬ 
quency  vircator.  It  optimizes  microwave  power  and 
efficiency  while  generating  high  frequency,  conerent 
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radiation.  For  low  frequency  operation  a  foil  diode 
or  reflex  triode  utilizing  a  high  transparency  mesh 
for  the  anode  can  be  used.  No  axial  magnetic  field 
should  be  employed  in  order  to  minimize  electron 
reflexing. 

A  schematic  of  the  configuration  being  used  In 
the  vircator  experiment1 7  at  Mission  Research  Corpor¬ 
ation  is  presented  in  Fig.  2.  The  pulse  power  para¬ 
meters  are  50  kv,  66  a  matched  impedance  and  a  pulse 
length  of  1  nsec.  The  magnetic  field  coil  is  capable 
of  attaining  60  kG.  We  anticipate  operating  at  fre¬ 
quencies  as  high  as  100  GHz.  Even  with  a  low  effi¬ 
ciency  of  2.51  we  will  produce  1  MW  of  RF  power.  A 
multi -channel  microwave  grating  spectrometer  in  the 
30-110  GHz  region  developed  by  MRC  will  be  the  chief 
diagnostic  on  this  experiment.17 

In  summary,  the  vircator  has  the  potential  for 
producing  very  high  power  microwave  pulses  in  tne 
centimeter  and  millimeter  wavelength  regimes.  In  a 
foil  less  diode  configuration  it  is  tunable  by  adjust¬ 
ing  the  imposed  axial  magnetic  field.  In  a  foil 
diode  or  reflex  triode  tuning  is  accomplished  by 
changing  the  A-K  gap  spacing.  The  microwave  genera¬ 
tion  will  be  coherent  and  efficient,  if  electron 
reflexing  into  the  diode  region  is  prevented. 

For  an  injected  current  v  >  3  v  ,  fQS£>“iD  //2*. 
Because  the  oscillating  beam  ts  equivalent  ro  a  de- 
fo’-mable  dipole,  the  preferred  waveguide  mode  for  an 
aosynmetric  beam  in  a  straight-walled  cylindrical 
guide  is  TM0n  where  n  =  D/x0,  D  is  the  waveguide 
diameter  and  x0  is  the  free  space  wavelength. 
Thus,  D/X  o  should  be  chosen  to  be  close  to  an 
integer  value.  Once  n  is  known,  the  phase  velocity, 
g'oup  velocity,  wavelength  and  impedance  of  the  wave 
in  the  guide  are  determined.  Similar  considerations 
hold  for  a  rectangular  waveguide.  Note  that  both 
vt  and  the  cutoff  wavelength,  xc,  depend  on  the 
guide  dimensions  and  geometry.  Both  must  be  consid¬ 
ered  in  choosing  an  experimental  configuration. 
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Abstract 

The  virtual  cathode  oscillator  (Vircator)  has  the 
potential  for  producing  very  high  power  Microwave  pul* 
ses  in  the  centimeter  and  millimeter  wavelength 
regimes.  In  a  foil  less  diode  configuration  it  is  tun¬ 
able  by  adjusting  the  Imposed  axial  magnetic  field. 
This  permits  high  frequency  operation.  In  a  foil 
diode  or  reflex  triode  tuning  is  accomplished  by 
changing  the  A-K  gap  spacing.  The  microwave  genera¬ 
tion  will  be  coherent  and  efficient,  if  electron 
reflexing  into  the  diode  region  is  prevented.  Giga¬ 
watt  power  levels  have  already  been  produced  in  the 
centimeter  regime. 

For  an  Injected  current  a  factor  of  three 
greater  than  the  limiting  current,  *osc  '"p*1  //?». 
Because  the  oscillating  beam  is  equivalent  to  a  defor¬ 
mable  dipole,  the  preferred  waveguide  mode  for  an 
axisymmetric  beam  in  a  straight -wal led  cylindrical 
guide  is  TMgn  where  n  •  D/a0,  D  Is  the  waveguide 
diameter  and  x0  is  the  free  space  wavelength.  Thus, 
0/xo  should  be  chosen  to  be  close  to  an  Integer 
value.  Note  that  both  the  limiting  current  and  the 
cutoff  wavelength  depend  on  the  guide  dimensions  and 
geometry.  Both  must  be  considered  in  choosing  an 
experimental  configuration. 

Introduction 

Of  the  several  millimeter  sources  that  are  in 
various  stages  of  development,  the  virtual  cathode 
oscillator  (Vircator)  has  a  combination  of  character¬ 
istics  which  recommend  it  for  high  frequency  use. 
First,  the  frequency  of  the  vircator  is  tunable  by 
changing  the  magnitude  of  an  Imposed  axial  magnetic 
field,  eliminating  any  requirement  to  change  the  phy¬ 
sical  structure  of  the  device.  A  single  vircator  will 
be  tunable  over  an  order  of  magnitude  in  frequency 
(e.g.,  10  GHz  -  100  GHz).  Second,  the  bandwidth  of 
the  generator  can  be  narrow  or  broad  based  on  magnetic 
field  shaping  and  the  use  of  beam  collimators  des¬ 
cribed  below.  Third,  because  the  vircator  functions 
above  the  space-charge  limiting  current  for  the  elec¬ 
tron  beam,  given  efficient  operation,  it  should  be 
capable  of  auch  higher  power  than  other  microwave 
sources.  Finally,  the  lack  of  passive  resonating 
structures  to  produce  the  transmitted  wave  reduces  the 
problem  of  field  emission.  This  also  Increases  the 
maximum  possible  generator  power. 


Experimentally,  the  virtual  cathode  has  already 
proven  Itself  to  be  a  copious  microwave  source  [1-7]. 
With  the  exception  of  the  Tomsk  experiment  [4]  the 
frequency  spectra  have  all  had  a  broad  bandwidth  and 
relatively  low  efficiency.  Nevertheless,  even  at  low 
efficiency  the  experiment  at  Marry  Diamond  Labora¬ 
tories  [5]  using  a  foilless  diode  (1  MV,  30  kA)  pro¬ 
duced  gigawatts  of  power  in  the  Ku  band.  It  is  one  of 
the  most  powerful  centimeter  wavelength  microwave 
sources  available.  It  will  be  shown  later  that  the 
foil  less  diode  in  a  shaped  axial  magnetic  field  repre¬ 
sents  the  best  configuration  for  a  high  frequency 
device. 

Microwave  Generation 

Although  only  scaling  relations  are  presently 
available,  qualitative  dependencies  of  virtual  cathode 
parameters  on  beam  kinetic  energy  and  Injected  current 
are  known  [8,9].  First,  potential  amplitude,  posi¬ 
tion,  and  oscillation  frequency  all  have  the  same 
functional  dependence  on  Injected  beam  current.  These 
parameters  asymptotically  approach  a  limiting  value 
for  current  above  the  space-charge  limit.  Second,  the 
fundamental  oscillation  frequency  is  approximately  the 
relativistic  beam  plasma  frequency  given  by 


where  n£  is  the  electron  beam  number  density  at  injec¬ 
tion  y0  is  the  beam  relativistic  factor,  e  is  the 
electron  charge  and  m  is  its  mass.  In  particular, 
the  oscillation  frequency  from  one-dimensional  elec¬ 
trostatic  and  two-dimensional  electromagnetic  numeri¬ 
cal  simulations  varies  such  that 

wpel  <  *  <  /7T  u.rel  •  (2) 

p  osc  ~  p 

The  value  of  /Zw  is  an  empirical  result  which  has  not 
yet  been  derived  theoretically.  The  value  of  »osc 
Increases  with  current  monotonically.  If  the  injec¬ 
tion  current  exceeds  the  space-charge  limiting  current 
by  a  factor  of  three  or  greater,  ugjc  is  close  to 
the  maximum  value.  Eq.  (1)  in  conjunction  with  the 
Child-Langmulr  law  describing  space-charge  limited 
diode  emission  explains  the  experimental  linear  depen¬ 
dence  of  frequency  on  the  square  root  of  diode  voltage 
in  foil  diodes  and  reflex  triodes.  This  analysis  also 


shows  that  the  freuqency  Is  Inversely  proportional  to 
the  anode-cathode  (A-K)  gap  spacing  0. 

The  virtual  cathode  oscillates  stably  at  a  set 
frequency  In  both  tine  and  space.  This  fluctuating 
potential  barrier  acts  as  a  gate  to  reflect  some  elec¬ 
trons  and  transmit  others.  The  motion  of  the  gate 
bunches  charge.  In  two  dimensions  the  charge  bunch 


probe  placed  between  the  anode  and  virtual  cathode. 

To  •  3.5,  u0  •  3.4  vt. 

By  analogy  this  configuration  represents  an  1C 
oscillator.  The  virtual  cathode  acts  as  a  capacitor 
to  store  the  beam  kinetic  energy.  During  that  portion 
of  the  limit  cycle  In  which  the  potential  Is  greater 
than  the  Injected  beam  energy,  charge  is  constrained 
to  remain  near  the  anode.  This  starves  the  virtual 
cathode  so  that  Its  amplitude  decreases  below  (y0 
-1)  mc2/e.  Once  this  occurs  the  charge  bunch  Is 
transmitted.  The  electron  motion  represents  a  large 
time  varying  current  through  an  Inductor.  The  pres¬ 
ence  of  charge  away  from  the  anode  reestablishes  the 
virtual  cathode,  and  the  cycle  repeats.  The  effect  on 
beam  current  can  be  examined  using  simulations.  A  net 
current  diagnostic  Is  given  in  Fig.  1  where  the  probe 
Is  positioned  between  the  anode  and  virtual  cathode. 
The  Injected  current  In  units  of  mcJ/e,  v0.  Is  3.4 
times  the  limiting  current,  vt,  and  y0  1$  3.5. 

Note  that  this  virtual  cathode  can  actually  reverse  the 
direction  of  current.  The  average  current  value  is 
uf 

The  oscillating  current  generates  microwaves 
[10,11],  The  wave  frequency  Is  the  oscillation 
frequency  of  the  virtual  cathode.  The  wave  propagates 
down  the  drift  tube  In  a  TM  waveguide  mode,  which 
determines  the  wavelength  and  phase  velocity  of  the 
wave  in  the  guide.  The  field  configuration  is  evident 


In  simulations  trfiere  there  Is  no  axial  magnetic  field. 
If  a  cold  beam  is  injected  and  azimuthal  symmetry  is 
assumed  by  the  code,  the  only  nonzero  fields  are  E,, 
Er,  and  Be  In  cylindrical  geometry.  These  three 
fields  define  a  TM  wave  traveling  In  the  z  direction. 

The  largest  Impediment  to  constructing  an  effi¬ 
cient  vlrcator  Is  the  effect  which  heating  has  on  the 
microwave  generation  efficiency.  As  noted  earlier 
[12],  beam  temperature  significantly  damps  out  the 
amplitude  of  the  potential  oscillation.  This  can  be 
understood  In  the  following  way.  For  a  monoenergetlc 
beam  all  of  the  charged  particles  bunch  at  the  same 
location.  Mathematically  this  represents  a  singu¬ 
larity  where  the  charge  density  goes  to  Infinity.  In 
reality  the  charge  bunch  Is  not  Infinitely  dense,  but 
It  does  become  several  times  greater  than  the  Initial 
beam  Injection  density.  The  severity  of  the  charge 
bunching  leads  to  efficient  microwave  generation.  If, 
on  the  other  hand,  the  beam  has  a  spread  In  axial 
momentum,  the  electrons  will  stop  at  different  loca¬ 
tions  In  the  potential  well.  This  tends  to  limit  the 
charge  bunching  and  the  amplitude  of  the  oscillating 
electric  and  magnetic  fields.  The  effect  of  beam 
temperature  In  reducing  the  RF  efficiency  of  the 
vlrcator  has  been  witnessed  In  one-dimensional  elec¬ 
tromagnetic  simulations  [13].  A  beam  spread  of  less 
than  3%  In  energy  reduces  the  microwave  generation 
efficiency  from  201  to  approximately  21.  Under  these 
conditions  the  vlrcator  Is  nothing  more  than  a 
Barkhausen  oscillator  [14].  Indeed,  the  low  effi¬ 
ciency  and  broad  bandwidth  observed  In  most  of  the 
experiments  to  date  can  probably  be  attributed  to  the 
effects  of  electron  reflexing  In  the  diode  region 
resulting  In  beam  heating. 

Vlrcator  Design 

Based  on  our  current  theoretical  knowledge  of  the 
virtual  cathode,  numerical  simulations,  and  experimen¬ 
tal  results,  the  following  points  must  be  considered 
In  the  design  of  a  coherent,  high  frequency  vlrcator. 
First,  generation  of  high  microwave  frequencies 
requires  large  beam  densities.  Extremely  high  beam 
densities  (n°>  lO1"  cm*J)  have  been  obtained  from  a 
foil  less  diode.  The  beam  plasma  frequency  scales  lin¬ 
early  with  the  electron  cyclotron  frequency  due  to  the 
Imposed  axial  magnetic  field  [15].  This  is  of  signi¬ 
ficance,  because  It  Implies  that  a  single  vlrcator  can 
be  tuned  over  an  order  of  magnitude  In  frequency  (e.g., 
10-100  GHz)  simply  by  changing  the  axial  magnetic 
field  strength  without  changing  the  physical  structure 
of  the  device. 
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Second,  both  the  oscillation  frequency  end  net 
current  asymptotical ly  approach  a  value  as  injected 
current  is  Increased  above  the  space-charge  limit. 
Thus,  using  a  very  large  value  of  v0/vt  does  not 
substantial ly  Increase  frequency  or  RF  efficiency 
where  efficiency  is  given  by 


eoa 


(y0-1)  me 


(3) 


A  foilless  diode  in  a  strong  axial  magnetic  field  pro¬ 
duces  a  very  thin  annular  beam.  Since  vt  for  an  an¬ 
nular  beam  is  larger  than  for  a  solid  beam  of  the  same 
area,  the  value  of  v0/vt  will  be  smaller  for  the 
same  beam  current. 


In  order  to  assure  narrow  bandwidth,  high  effi¬ 
ciency  microwave  generation  at  high  or  low  frequency 
from  the  vircator,  the  following  characteristics  must 
be  met.  First,  no  reflexing  of  electrons  in  the 
region  between  the  anode  and  real  cathode  should 
occur.  Any  axial  magnetic  field  must  be  shaped  to 
divert  the  electrons,  or  flux  excluders  must  be 
employed  to  confine  the  magnetic  field  to  the  diode 
region.  The  latter  arrangement  will  allow  the  radial 
space-charge  electric  field  to  perform  the  role  of 
expelling  electrons  to  the  waveguide  wall.  In  addi¬ 
tion,  if  the  beam  is  annular,  a  collimator  may  be  used 
to  help  prevent  reflexing  of  electrons  back  to  the 
cathode.  Second,  the  electron  beam  nust  be  cold. 
Experimental  [16]  and  theoretical  [15]  results  indi¬ 
cate  that  foil  less  diodes  create  low  emittance  beams. 
Laminar  flow,  where  the  electron  Larmor  orbit  Is 
smaller  than  the  beam  thickness,  is  obtained  when  [16] 

.1/2  c 


-c  >  <*o  * 


/  a«/2 


(«) 


where  a  is  the  orbit  radius,  6  is  the  radial  spacing 
between  the  cathode  and  drift  tube  wall  (which  acts  as 
the  anode)  and  uc  is  the  electron  cyclotron  fre¬ 
quency  given  by  eBz/mc.  Low  beam  scatter  is  also 
assured,  because  of  the  lack  of  a  foil.  Finally,  the 
diode  voltage  and  Injected  current  must  be  constant, 
•tore  appropriately  stated,  the  Impedance  must  be  con¬ 
stant.  Flat-top  voltage  pulses  can  be  attained  in  a 
variety  of  ways  in  several  diode  configurations.  How¬ 
ever,  at  high  voltages  the  foil  less  diode  operates  as 


a  purely  resistive  load,  therefore,  «osc*  /T7?  *  1//7 


is  constant.  Also,  absence  of  diode  closure  in  some 
foilless  diode  experiments  makes  a  long  pulse  device 
possible. 


It  is  evident  from  this  discussion  that  the  foi 1  - 
less  diode  In  a  strong  axial  magnetic  field  represents 
the  optimal  configuration  for  a  high  frequency 


vircator.  It  optimizes  microwave  power  and  efficiency 
while  generating  high  frequency,  coherent  radiation. 
For  low  frequency  operation  a  foil  diode  or  reflex 
triode  utilizing  a  high  transparency  mesh  for  the 
anode  Is  the  optimal  design.  No  axial  magnetic  field 
should  be  employed  in  this  case  in  order  to  minimize 
electron  reflexing. 

Conclusion 

In  suenary,  the  vircator  has  the  potential  for 
producing  very  high  power  microwave  pulses  in  the  cen¬ 
timeter  and  millimeter  wavelength  regimes.  In  a  foi  1  - 
less  diode  configuration  It  is  tunable  by  adjusting 
the  Imposed  axial  magnetic  field.  In  a  foil  diode  or 
reflex  triode  tuning  is  accomplished  by  changing  the 
A-K  gap  spacing.  The  microwave  generation  will  be 
coherent  and  efficient,  if  electron  reflexing  into  the 
diode  region  is  prevented. 

For  an  injected  current  vQ>  3  v_,  the  frequency 
is  approximately 

f  ■  10.2Vj(kA/cm2)/6Yo  GHz  (5) 

where  t  is  v/c.  Because  the  oscillating  beam  is  equi¬ 
valent  to  a  deformable  dipole,  the  preferred  waveguide 
mode  for  an  axisymmetric  beam  in  a  stralght-wal led 
cylindrical  guide  is  TMgp  where  n  *  D/x0,  D  is  the 
waveguide  diameter  and  x0  is  the  free  space  wave¬ 
length.  Thus,  D/x0  should  be  chosen  to  be  close  to 
an  Integer  value,  in  order  to  be  near  cutoff  for  that 
waveguide  mode.  This  maximizes  coupling  to  the  no-i- 
resonant  waveguide,  because  E,/B  which  is  equal  to 
the  wave  phase  velocity  is  substantially  larger  than 
c.  Once  n  Is  known,  the  phase  velocity,  group  velo¬ 
city,  wavelength  and  impedance  of  the  wave  in  the 
guide  are  determined.  Similar  considerations  hold  for 
a  rectangular  waveguide.  Note  that  both  vz  and  the 
cutoff  wavelength,  kc,  depend  on  the  guide  dimen¬ 
sions  and  geometry.  Both  must  be  considered  in 
choosing  an  experimental  configuration. 
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Effects  of  thermal  spread  on  the  space  charge  limit 
of  an  electron  beam 
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An  asymptotic  analysis  is  carried  out  to  calculate  the  effects  of  a  small  thermal 
spread  in  the  injection  energy  of  an  electron  beam  on  its  space  charge  limit.  It 
is  found  that  the  space  charge  limit  is  lowered  proportionally  to  the  beam  tem¬ 
perature  T  near  T  =  0. 


1.  Introduction 

Recent  applications  of  intense  charged  particle  beams  in  such  areas  as  inertial 
confinement  fusion  and  microwave  generation  (Coutsias  &  Sullivan  1983.  and 
references  therein)  has  necessitated  a  deepening  of  our  understanding  of  the 
basic  physics  of  space  charge  limited  flows.  An  excellent  review  of  our  present 
state  of  knowledge  can  be  found  in  Miller  (1982). 

In  particular,  many  authors  have  carried  out  calculations  of  the  space  charge 
limit  (SCL)  of  electron  and  ion  beams  (Voronin.  Zozulya  &  Lebedev  1972;  Read 
&  Nation  1975;  Genoni  &  Proctor  1980)  These  are  concerned  mostly  with  mono- 
energetic  beams  in  various  geometries.  Although  the  importance  of  thermal 
effects  is  recognized,  no  analytical  estimates  of  the  effect  of  thermal  spread  on 
the  beam  kinetic  energy  at  injection  have  appeared. 

Here  we  present  an  asymptotic  method  to  estimate  the  modification  of  the 
SCL  due  to  a  small  thermal  spread.  For  simplicity  we  limit  our  discussion  to 
one-dimensional,  classical  motion.  However,  the  method  can  be  applied  to  any 
of  the  other  situations  for  w  hich  SCL  estimates  exist  and  produce  appropriate 
corrections. 


2.  Effect  of  temperature  on  the  space  charge  limit 

An  adequate  discussion  of  the  properties  of  a  non-relativistic  electron  beam  in 
one  dimension  can  be  found  in  Coutsias  A  Sullivan  (1983).  Here  we  shall  treat 
the  beam  as  a  one-dimensional  electron  gas  flowing  between  two  conducting 
grid  planes,  at  fixed  potential.  For  the  density  range  we  consider  (i  1  O’3  cm-*) 
the  usual  collisionless  approximation  is  valid,  and  thus  the  electron  distribution 
function  satisfies  the  Vlasov  equation 

dJ+u?xf+(e/m)EdJ=0  (1) 
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where  t,  x,  u  are  the  time,  position  and  velocity  variables,  respectively,  and  E 
is  the  self-consistent  electric  field. 

The  electron  charge  density  is  given  by 


n  =  J  f(x,u,t 


and  the  electric  field  is  found  from 

dxE  =  (l/e0)n,  (3) 

while  the  potential  $(x,t)  is  given  by 

dx$  =  -E.  (4) 

The  boundary  conditions  are  given  for  the  potential 

<5(0,f)  =  <50>0,  o.  (5) 

and  the  distribution  function 

f(0,u,t)  specified  for  u>0  (incoming  flow  at  x  =  0)  l 
/(/,  a,0  =  0  for  u  <  0  (no  incoming  flow  at  x  =  /)/ 

This  specification  is  valid  also  in  the  presence  of  multiple  streams. 

We  shall  model  the  effect  of  finite  source  temperature  by  specifying  the  in¬ 
coming  distribution  f(0,u,t),  a  >  0  as 


where  T,  assumed  to  be  small,  plays  the  role  of  an  effective  beam  temperature 
For  small  enough  T.  we  shall  assume  that  we  have  a  regime  of  steady  behaviour, 
in  analogy  to  the  cold  beam  case.  For  steady  states,  particle  energy  is  conserved 
and  the  solution  to  (1)  is  of  the  form  (Davidson  1974) 

m[V-  («*  +  (2e/r7i)  (<5(r)  —  <50))*1 


Then,  combining  (2).  (3),  (4)  and  (8)  we  find  that  the  potential  satisfies  the 
equation 

-  »■  <•> 

where  um,  the  velocity  cut-off.  is  equal  to 

um(*)  =  ±[(2e/wi)(<Jm-<5(or))]*.  (10) 

Here  is  the  potential  minimum,  the  ( +  )  sign  applies  to  the  right  and  the  (  -  ) 
to  the  left  of  the  position  j  =  £  of  the  potential  minimum  as  shown  in  Appendix  A. 
We  define 

So  (9)  <an  be  written  as 


,2; 
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By  introducing  the  variable 

e  =  (u*+(2e/m)($5(x)-$J0))*- V  (13) 

the  integral  ;  T)  in  (12)  can  be  rewritten  as 

The  first  of  the  integrals  in  (14)  represents  the  transmitted  flow,  while  the  Becond 
is  due  to  particles  without  sufficient  energy  to  cross  the  potential  minimum  and 
which  are  therefore  reflected  and  return  to  the  anode.  Equation  (12)  with  the 
integral  term  given  by  (14)  is  very  hard  to  solve  for  a  finite  temperature  T,  but 
for  small  T  we  can  approximate  the  integrals  in  (14)  by  Laplace’s  method.  As 
is  well  known  (Erdelyi  1956),  in  approximating  integrals  of  this  type  with  a 
strong  maximum  at  an  interior  point,  the  dominant  contribution  comes  from 
the  neighbourhood  of  this  point.  Thus,  in  (14)  we  can  approximate  7(0;  T)  by 

/(f>; T>  -  eip  ( ■ -  m ) iu + 0  (e,p  ( ■ -  irl))  1 ,4“, 

where  —sm  =  l’- ((2e/m)  (<?>m  -  $50))l.  Then  provided  V  is  large  enough  so  that 
the  mean  energy  of  the  beam  is  never  of  order  0(kT).  the  correction  term  goes  to 
zero  faster  than  any  power  of  T  as  T  -*■  0,  and  therefore  is  negligible  to  the  order 
that  we  carry  the  calculations.  Nevertheless,  it  gives  us  an  estimate  of  the  domain 
of  validity  of  the  subsequent  discussion,  for  which  we  need 

exj>(-msl,/2kT)  <g  kT.  (146) 

Therefore  we  write 

f*  /  msi\,  f®  (s+  V)e*p(-msi/2kT) 

..‘M-zrr-J <15’ 

and  we  get  for  the  potential  the  approximate  equation 

A  ,  "of  m  V  f"  (s  4-  Tr) exp ( - ms2/2kT) 

<*"  +  e0Wr)  J_.  ((s+r)*-(2e/ni)(^)-?i0))irf  (  ) 

This  can  be  integrated  onoe  to  give 

+  (IT) 


where 


x  ^(s+  F)1- ((*(*)-&>))  exp ( —  m3t/2kT)ds.  (18) 


Integrating  once  more  we  find 


=  S ±S(<(>m,<f>;T ). 


Imposing  the  boundary  conditions  at  x  =  0  and  x  =  l  we  are  led  to  the  two 
equations  determining  £  and  <f>m: 

0  l  =  £  +  S($m,0;T).  (20) 
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Eliminating  £  we  find 

l  =  Stfm.tu,T)  +  S(tm,0-,T)  (21) 

which  must  be  analysed  in  order  to  determine,  among  other  things,  the  desired 
steady  state  for  the  potential  $5  (and  hence  the  distribution  function/)  and  the 
SCL  for  small  temperature  T  >  0. 

To  get  an  approximate  expression  for  (j>m  as  T  ->  0  we  note  that  J($\ T)  can 
be  approximated  by  Laplace's  method  if  we  expand  the  non-exponential  part 
of  the  integrand  in  a  Taylor  series  about  s  =  0  and  integrate  term  by  term.  We 
find,  after  some  algebra  that 


I 


J($;T)±  VR($)  + 


eV  (I?2($)-(e/m)(<t>-<f>0) 


Tt 


R 3 


+  0(e2) 


(22) 


where  e  =  (2kT/m)4,\  and  R(<f>)  =  (Vi-(2e/m)(<f>(x)-(f)0))i.  Using  (22)  we 
now  approximate  5 : 


=  2-i  P 
k 


d<t> 


C.qf.-V 

\‘2n0m  1  / 


\-«0' 

d<i>  (,  e  l(r2-<t>)Rl-(R2m 


(R-RJ* 


(‘-if1 


R’RUR-RJ 


J+0(e2))  (23) 


where  we  set  R  =  R($).  Rm  =  R(<t>m).  <i>  =  (f/ro)(0(x)-0o),  <I>m  =  (e/m)($m-$u) 
and  since  d$  =  -  (m/e)RdR. 


sit*.-* T)  ^  m) -  *(&,,»*  m) - 2 Rtfm))  +£(2t^)1 


((/?  Rj  (4^ 


r* 

4  Ri 


3T2 

8  Ri, 


Ri 

Ri 


T*  \  3l’2  Ri\ 

md-^-'iz)+0^  1241 


This  can  be  substituted  in  (21 )  to  find  (f>m  which  in  turn  will  allow  us  to  determine 
£,  the  position  of  the  potential  minimum  from  (20)  and.  finally,  we  can  combine 
all  this  information  in  (19)  to  get  the  desired  relation  between  x  and 

For  simplicity  we  demonstrate  this  for  the  unbiased  case,  <f>a  =  0.  Letting 


and  introducing 
(21 )  becomes 


i?(0)  =  Rtf0)  =  1\  s  =  RJV, 
a  =  9cn0?2/8£0mr?, 


ai 


(1+2  s)  + 


3  _J _ 1_ 

45  45*  85s 


e  =  e/n, 


3  sec-1  s-* 
8  si  ( 1  —  s)i 


(25) 


We  note  that  to  leading  order  (c  =  0).  (24)  reduces  to  the  usual  expression  for 
cold  beams.  By  including  the  0{r )  corrections  we  can  find  the  first-order  correc¬ 
tion  to  the  SCL.  This  is  the  value  of  a  for  which  (24)  has  a  double  root.  We  find 


a  =  o0  +  eo,  +  s  =  5„  +  £S,  +  ..., 
and  substitution  yields 


au  —  2,  s0  —  J, 
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Figure  1.  Phase  plane,  (a)  Transmitted  pamoles,  (b)  reflected  particles,  (c)  um  =  ((2 e/m) 
(<j>m  —  <p\ (</)  um  ~  —  ((2c/m)  —  0(.r)))i,  (e)  £,  location  of  potential  minimum, 


front  which  wc  find  o,: 

Oj  =  —  3(1  +  Jff)  -  —  10  07 
That  is,  the  SCL  is  given  by 


enj 2 


so  that,  for  T  4  1,  the  SCL  is  decreased  in  proportion  to  the  temperature. 


This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research  under 
contract  no.  AFOSR-82  0277  and  F49C20-82-C-0014 


Appendix  A 

From  the  steady-state  momentum  equation 

uux=  (e/m)E  (A  1 ) 

we  find  by  integration  that 

lu'  +  —p(r)  =  constant.  (A  2) 

m 

At  x  =  0 ,  u  =  u0.  p  =  p0  and 

*M2  +  mp(‘r)  =  +  (A  3) 

Since  (1)  is  a  statement  of  the  conservation  of  /  along  particle  trajectories,  we 
find  that  if  u  and  u0  are  related  by  (A  3)  then 

f{r,u.t)  =/(0,  u0.t)  (A  4) 

and  expression  (8)  follows. 

A  complication  arises  when  we  are  dealing  with  reflected  particles  since  u  in 
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(A  3)  can  have  either  sign.  Assuming  that  the  potential  achieves  a  unique  mini¬ 
mum  p  =  pm  at  some  location  x  =  £,  we  find  that  particles  that  arrive  at  £  with 
zero  velocity  must  start  at  x  =  0  with  velocity  um(0)  such  that 

i“m(°)  +  ^o  =  + 

At  each  point  x  <  £  particles  with  velocities  less  than  vm  will  be  reflected  before 
reaching  x  =  £  while  faster  particles  will  be  transmitted  across  the  potential 
barrier.  To  the  left  of  the  potential  minimum  therefore,  we  shall  have  particles 
with  velocities  larger  than  ((2c/m)  (pm-  p(x)))l  composing  the  flow  that  will  be 
transmitted,  and  particles  with  velocities  in  the  range 

|«|  <  (( 2e/m)(pm-p(x)))l 

composing  the  counterstreaming  flow  (figure  1).  To  the  right  of  the  potential 
minimum  we  have  only  particles  with  velocities  larger  than  um,  so  that  the 
potential  minimum  has  filtered  out  of  the  flow  particles  whose  energies  were  too 
small  to  traverse  it. 
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